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I. INTRODUCTION 

There has been a vigorous and global program of ex- 
periments and theoretical developments in the last 25 
years aimed at understanding the internal spin structure 
of the proton. How is the proton's spin built up from 
the spin and orbital angular momentum of the quarks 
and gluons inside? Tremendous progress has been made 
with unraveling the proton's spin structure with advances 
in experimental techniques, theoretical models, pertur- 
bative QCD, non-perturbative QCD and lattice calcula- 
tions. 

This activity was inspired by the initial European 
Muon Collaboration (EMC) data which suggested the 
puzzling result that quark intrinsic spin co ntributes little 
to the proton's spin (jAshman et al , 19881 ). Today there 
is good convergence of the theoretical and experimen- 
tal understanding the proton's longitudinal spin struc- 
ture. Further puzzling data in measurements of trans- 
verse single-spin asymmetries revealed up to 40% asym- 
metries in proton-proton collisions (and 5-10% in lepton- 
nucleon collisions with unpolarized leptons and trans- 
versely polarized nucleons) which persist to high energies. 
These single-spin asymmetries indicate significant spin- 
orbit coupling in the nucleon associated with quark trans- 
verse momentum and the bound state structure of the nu- 
cleon. The study of transverse momentum and associated 
orbital angular-momentum processes has spawned new 
programs to map out the three-dimensional structure of 
the nucleon. In this article we review these developments 
highlighting the considerable and exciting developments 
in QCD spin physics in recent years, together with an 
outlook to the future: What are the main open questions 
and the planned experiments to help answer them? 

In 1988 EMC published their polarized deep inelastic 
measurement of the proton's gi spin dependent structure 
function and the flavor-singlet axial-charge g^ (the nu- 
cleon's "quark spin content" ) suggesting that quark spins 
summed over up, down and strange quark flavors con- 
tribute only a small fraction of the proton's spin. This 
result inspired considerable theoretical activity and new 
experiments at CERN, SLAC, DESY, Jefferson Labo- 
ratory (JLab) and the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory (BNL) to 
understand the spin structure of the nucleon. The first 
task was to check the initial curious result from EMC and 
second to resolve the spin-flavor structure of the proton. 
How is the spin content of the proton distributed among 
the valence and sea quarks and gluons? What about or- 
bital angular momentum in the nucleon? 

We now know that the nucleon's flavor-singlet axial- 
charge measured in polarized deep inelastic scattering is 
9a ~ 0.35. This value was surprising from the view- 
point of early quark models. In the static quark model 
- the eightfold way picture of Cell-Mann - before in- 



clusion of quark motion, quark spin contributes 100% 
of the proton's spin. Relativistic quark models without 
gluonic or pion cloud degrees of freedom generally pre- 
dict about 60% of the proton's spin should be carried 
by the quarks, with the remaining 40% in quark orbital 
angular momentum. Today data and theory point to a 
consistent picture where the proton spin puzzle is a va- 
lence quark effect. Valence quark contributions to g^ 
approximately saturate the measured value. While po- 
larized glue may contribute a significant fraction of the 
proton's spin (perhaps up to 50% at the scale of present 
experiments), sea quark and QCD gluon corrections to 
the singlet axial-charge are small and within the expecta- 
tions of quark models. The pion cloud of the nucleon acts 
to shift angular momentum from spin to orbital angular 
momentum and induces SU(3) breaking in the nucleon's 
axial-charges. There is also a fascinating theoretical pos- 
sibility that the valence quarks may polarize the QCD 
vacuum in a nucleon through gluon topological effects so 
that some fraction of the proton's singlet axial-charge re- 
sides at zero parton momentum (or Bjorken x). Non-zero 
orbital angular momentum of the valence quarks is ex- 
pected, induced also by confinement which introduces a 
transverse scale in the physics. This orbital angular mo- 
mentum through spin-orbit coupling is a prime candidate 
to explain the large single spin asymmetries observed in 
proton-proton collisions. Information about quark total 
angular momentum in the proton can be extracted from 
deeply virtual Compton scattering and high-energy sin- 
gle spin asymmetry data in model-dependent analyses. 
The results are consistent with QCD lattice calculations. 

This Review is organized as follows. In the first part 
(Sections II— III) we give a brief introduction to nucleon 
spin physics and the experiments that have been per- 
formed to investigate it. Then, in Section IV, we discuss 
the proton spin puzzle and the small value of g^ ex- 
tracted from polarized deep inelastic scattering. In Sec- 
tion V we give an overview of the present global pro- 
gram aimed at disentangling the spin-flavor structure of 
the proton. Section VI covers the theoretical interpreta- 
tion of longitudinal spin data and understanding of the 
proton spin puzzle. We next turn our attention to the 
transverse structure of the nucleon and manifestations 
of orbital angular momentum in the nucleon in Section 
VII. This discussion introduces generalized parton dis- 
tributions (GPDs), which describe hard exclusive reac- 
tion processes, and transverse momentum dependent dis- 
tributions (TMDs), which describe spin-momentum cor- 
relations and spin-orbit couplings in the nucleon. The 
TMDs are manifest in high-energy single-spin and az- 
imuthal asymmetries. A summary of key issues and chal- 
lenging questions for the next generation of experiments 
is given in Section VIII and IX. 

Earlier review articles on the spin structure of 
the proton as well as complementary more recent 
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II. SPIN STRUCTURE FUNCTIONS AND PARTON 
DISTRIBUTIONS 

Our knowledge about the high-energy spin structure 
of the nucleon comes from both polarized deep inelastic 
scattering (DIS) experiments and high-energy polarized 
proton-proton collisions. Polarized deep inelastic scat- 
tering (pDIS) experiments involve scattering a longitu- 
dinally polarized high-energy lepton beam from a longi- 
tudinally or transversely polarized nucleon at large mo- 
mentum transfer. Inclusive measurements, where only 
the scattered lepton is observed in the final state, and 
semi-inclusive measurements, where one tags on at least 
one high-energy final state hadron in coincidence with the 
scattered lepton, have been performed. The experiments 
were performed with an electron beam at SLAC and 
JLAB, with electron and positron beams at DESY and 
with muon beams at CERN. In proton-proton scattering 
the protons are either longitudinally or transversely po- 
larized. Polarized deep inelastic scattering experiments 
have so far all been performed using a fixed target. A fu- 
ture polarized electron-ion collider is in planning. Details 
of the experiments are given in Section III. Historically, 
information about the proton's internal spin structure 
came first from measuring the proton's g\ and gi spin 
structure functions in inclusive deep inelastic scattering 
and, more recently, from semi-inclusive reactions in both 
lepton-nucleon and proton-proton collisions and hard ex- 
clusive processes in lepton-nucleon scattering. 

Measurements with longitudinally polarized targets 
and beams tell us about the helicity distributions of 
quarks and gluons in the nucleon, which at leading or- 
der can be thought of as the difference in probability of 
finding a parton with longitudinal polarization parallel 
or anti-parallel to that of the nucleon. Measurements 
with transversely polarized targets are particularly sen- 
sitive to quark and gluon transverse and orbital angular 
momentum. Studies of transverse degrees of freedom in 
the nucleon and in fragmentation processes are a current 
subject of experimental investigation with sensitivity to 
spin-orbit couplings in QCD. 

For polarized lepton-proton scattering, specialize to 
the target rest frame and let E denote the energy of the 
incident lepton which is scattered through an angle 6 to 
emerge in the final state with energy E' . Let y! de- 
note the longitudinal polarization of the lepton beam. In 
photon-nucleon scattering the spin dependent structure 



functions g\ and $2 are defined through the imaginary 
part of the forward Compton scattering amplitude. The 
structure functions contain all of the target-dependent 
information in the deep inelastic process. Consider the 
amplitude for forward scattering of a photon carrying 
momentum (q 2 = —Q 2 < 0) from a polarized nucleon 
with momentum p M , mass M and spin s M . We work with 
the kinematic Bjorken variable x = Q 2 /2p ■ q = Q 2 /2Mv 
where v = p-q/M = E— E', and let y = p-q/p-k = v/E. 
For a longitudinally polarized proton target (with spin 
denoted ffJJ.) the unpolarized and polarized differential 
cross-sections are 



dxdy 
2ira 2 
MEx 2 y 2 
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where the mass of the lepton is neglected. The relation 
between the structure functions in deep inelastic lepton- 
nucleon scattering and the virtual-photon nucleon cross- 
sect ions is d i scusse d and derived in various textbooks, 
e.g. iRobertsI (|l990l ). One finds 



Ai = 



9i 



Q 



92 



El 

Fx 



(3) 



where a 3 and a 1 are the cross-sections for the absorption 
of a transversely polarized photon with spin polarized 
parallel and anti-parallel to the spin of the longitudinally 
polarized nucleon. For a longitudinal polarized target the 
gi contribution to the differential cross-section and the 
longitudinal spin asymmetry is suppressed relative to the 
gi contribution by the kinematic factor M/E -C I. For 
a transverse polarized target this kinematic suppression 
factor for g2 is missing implying that tra nsverse polariza - 
tion is vital to measure 02- We refer to iRoberts (|l990h 
and Windmoldersl ( 2002 ) for the procedure how the spin 
dependent structure functions are extracted from the 
spin asymmetries measured in polarized deep inelastic 
scattering. 

In high- Q 2 deep inelastic scattering the structure func- 
tions F\, F2, g\ and gi exhibit approximate scaling. They 
are to a very good approximation independent of Q 2 and 
depend only on Bjorken x. (The small Q 2 dependence 
which is present in these structure functions is logarith- 
mic and determined by perturbative QCD evolution.) 
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In the (pre-QCD) parton model the deep inelastic 
structure functions F\ and F% are written as 

F 1 (x) = ±.F a {x) =\j24{Q + fiW W 



and the polarized structure function g\ is 
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(5) 



Here e q denotes the electric charge of the struck quark 
and 



{q + q}(x) 
Aq{x) 



(9 f 



f)(x) 



(q l 



t)(x) 
t)(x) 



(6) 



denote the spin- independent (unpolarized) and spin- 
dependent quark parton distributions which measure the 
distribution of quark momentum and spin in the proton. 
For example, cf(x) is interpreted as the probability to 
find an anti-quark of flavor q with plus component of 
momentum xp + (p + — po + P3 is the plus component 
of the target proton's momentum) and spin polarized 
in the same direction as the spin of the target proton. 
When we integrate out the momentum fraction x the 
quantity Aq = J dx Aq{x) is interpreted as the frac- 
tion of the proton's spin which is carried by quarks (and 
anti-quarks) of flavor q. Hence summing over the up, 
down and strange quark Aq contributions gives the total 
fraction of the proton's spin carried by the spins of these 
quarks. 

What values should we expect for the Aql First, con- 
sider the static quark model. The simple SU(6) proton 
wavefunction 

\P t) = -j=\u t Ms=o) + ^==|u t Ms=i) 



-~|u I (ud) s =i) - ~|d t (vm) s =i) 



(7) 



yields the values Au — Ad = | and Au + Ad — 1 . In rel- 
ativistic quark models one has to take into account the 



four-component Dirac spinor ip 



f 

ia-fg 



The lower 



component of the Dirac spinor is p-wave with intrinsic 
spin primarily poi nting in the opposite direc tion to the 
spin of the proton ( Jaffe and Manoharlll990() . Relativis- 
tic effects renormalize the axial charges by the depolar- 
ization factor 0.65 with a net transfer of angular momen- 
tum from intrinsic spin to orbital angular momentum. In 
QCD and in more sophisticated models further depolar- 
ization is induced by gluonic and pion-cloud degrees of 
freedom - see Section VI. 



In QCD the flavor-singlet combination of the 
Aq(x) quark parton distributions mixes with the 
spin-dependent gluon dist ribution under Q 2 evolution 
(jAltarelli and Parisii Il977l ) . This spin dependent gluon 
distribution measures the momentum and spin depen- 
dence of glue in the proton. The second spin structure 
function 172 vanishes without the effect of quark trans- 
verse momentum and has a non-t rivial parton interpre- 
tation (|jaffel . I199CA iRobertsl . Il99dh . 

The parton model description of polarized deep inelas- 
tic scattering involves writing the deep inelastic structure 
functions as the sum over the convolution of "soft" quark 
and gluon parton distributions with "hard" photon- 
parton scattering coefficients 



&{x) = J^(Au-Ad) 



-I (Au + Ad 



1 

36 



(Au + Ad - 2As) 



As)®C q s + fAg{z^ 2 )®C 9 



(8) 



Here Aq{x) and Ag(x) denote the polarized quark and 
gluon parton distributions, C q and C 9 denote the corre- 
sponding hard-scattering coefficients, and / is the num- 
ber of quark flavors liberated into the final state (/ = 3 
below the charm production threshold) . The parton dis- 
tributions contain all the target-dependent information 
and describe a flux of quark and gluon partons into the 
(target independent) interaction between the hard pho- 
ton and the parton which is described by the coefficients 
C q and C 9 . These coefficients are calculated using per- 
turbative QCD via the cross-section for the hard photon 
scattering from a quark or gluon parton "target" . They 
are independent of infra-red mass singularities (terms in- 
volving the quark mass or virtuality of the parton in the 
photon-parton collision) which are absorbed into the par- 
ton distributions (and softened by confinement related 
physics). If the same recipe ("factorization scheme") for 
separating hard and soft parts of the parton phase space 
is applied consistently to all hard processes then the fac- 
torization theorem asserts that the parton distributions 
that one extracts from different experiments are process 
independent. In other words, the same polarized quark 
and gluon distributions should be obtained from experi- 
ments involving polarized hard QCD processes in polar- 
ized proton-proton collisions and polarized deep inelastic 
scattering experiments. For example, colliding longitudi- 
nally polarized proton beams provides sensitivity to the 
gluon- helicity distribution function at leading order. For 
hadron production with transverse momentum px, the 
helicity-dependent difference in hadron production is de- 
fined as 



dAa 



da++ da+~ 



dpi 



(9) 



4 



where the superscripts ++ and H — refer to same and 
opposite helicity combinations of the colliding protons. 
Factorization allows this to be written as a convolution 
of the long- and short-distance terms summed over all 
possible flavors for the partonic interaction a + b — > jet + 
X 



dAa 

dpT 



J dx a dx b Af a (x a ,n)Af b (x b ,Li) 

ab 

dAa ab^jet+X 
X {X a Pa,X b P b ,fi). 

dpr 



(10) 



Here P a and P b denote the momenta of the incident pro- 
tons; f(x, /i) are the polarized parton distributions of the 
colliding partons carrying light-cone momentum fraction 
x evaluated at factorization and renormalization scale /i. 
The helicity-dependent difference in the cross-section of 
the hard partonic scattering a + b — > jet + X is denoted 
by dAa and is calculable in perturbative QCD. Partonic 
cross-section calculations are carried out to finite order 
in a s and have a dependence on factorization and renor- 
malization scales, denoted \x. The final hadronic cross- 
section is independent of the factorization and renormal- 
ization scales and the scheme used. The QCD parton 
model treatment readily generalizes to the production 
of high-energy hadrons in the final state, with the pro- 
duced "fast" hadron carrying a significant fraction of the 
momentum of a "parent" parton. The parton-to-hadron 
process is parametrized by fragmentation functions which 
also obey process-independent factorization in perturba- 
tive QCD calculations. 

Analogous to the helicity distributions measured with 
longitudinal polarization, transversity distributions de- 
scribe the density of transversely polarized quarks inside 



a tran sversely polarized proton, see e.g. iBarone et al. 
(2002). The transversity distributions, which were in- 



troduc ed i nlRalston and Soperl (|l979t ). Artru and Mekhfi 
(|l990h and l.Taffe and .Til (jl992h . are interpreted in parton 
language as follows. Consider a nucleon moving with 
(infinite) momentum in the {^-direction, but polarized 
transverse to £3. Then Sq(x) (also denoted Axq{x) and 
h\{x) in the literature) counts the quarks with flavor q, 
momentum fraction x and their spin parallel to the spin 
of a nucleon minus the number anti-parallel. That is, 
in analogy with Eq. ([5]), Sq(x) measures the distribution 
of partons with transverse polarization in a transversely 
polarized nucleon, viz. 

Sq(x) = q^x) + q^(x) - q l (x) - q l (x). (11) 

In a helicity basis transversity corresponds to hel icity-flip 
making it a probe of chiral symmetry breaking ( Collins! . 
1993(). There is no gluon analogue of transversity in 
the nucleon so Sq evolves in Q 2 like a valence or non- 
singlet quark distribution, without mixing with glue. 
If quarks moved non-relativistically in the nucleon Sq 



and Aq would be identical since rotations and Euclidean 
boosts commute and a series of boosts and rotations can 
convert a longitudinally polarized nucleon into a trans- 
versely polarized nucleon at infinite momentum. The 
difference between the transversity and helicity distribu- 
tions reflects the relativistic character of quark motion in 
the nucleon. 

Following the discovery that the quark spin contribu- 
tion to the proton's spin is small, there has been a vig- 
orous program to measure the separate contributions of 
up, down and strange quark flavors as well as the gluon 
spin and the orbital contributions. This has inspired 
dedicated spin programs in semi-inclusive deep inelas- 
tic scattering (SIDIS) and polarized proton-proton col- 
lisions to measure the separate valence and sea quark 
as well as gluon polarization. As efforts to investigate 
nucleon spin in more detail intensified and various ex- 
perimental programs were being developed in the 1990s, 
new theoretical ideas arose as well. TMD distributions, 
describing spin-momen tum correlati ons in the nucleon, 
were initially proposed (|SiversL ll990'l to explain the very 
large transverse single spin asymmetries involved in po- 
larized h adronic scatt e ring t hat w ere first observed in th e 
1970s by iKlem et all (Il976l) andlDraeoset efdl (Il978l) 

~ .17)11 <ll997bh 



The GPDs introduced in Mueller et al. (1994 



and iRadvushkin ( 1997f ) to describe hard exclusive reac- 
tions provided for the first time a means of describing 
the radial position distributions of partons at a specific 
longitudinal momentum within the nucleon. Both TMD 
distributions and GPDs offer links to the orbital angu- 
lar momentum contributions to the nucleon's spin. These 
processes and the present status of experimental and the- 
oretical investigation are described in Section VII. 



III. EXPERIMENTS 

Experiments that have probed the nucleon spin struc- 
ture are outlined in Table I. This includes both polar- 
ized deep inelastic lepton-nucleon scattering and proton- 
proton collision experiments. Considerable effort was in- 
vested in developing polarized beam and target technol- 
ogy, yielding physics results with ever increasing preci- 
sion. The first experiments focused on inclusive deep 
inelastic measurements of nucleon spin structure. More 
recent experiments, described in detail below, were able 
to detect and identify hadrons in the final state lead- 
ing to new probes of the nucleon in semi-inclusive and 
hard exclusive reactions. Future experimental programs 
(COMPASS-II, the 12 GeV upgrade of JLab and experi- 
ments at Fermilab and RHIC) with high luminosity and 
acceptance are planned to explore the three-dimensional 
structure of the nucleon in spatial and transverse mo- 
mentum degrees of freedom. We discuss these future 
programs in Section VIII. 
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TABLE I High energy spin experiments: the kinematic ranges in x and Q 2 correspond to the average kinematic values of the 
highest statistics measurement of each experiment, which is typically the inclusive spin asymmetry; x denotes Bjorken x unless 
specified. 



Experiment 


Year 


Beam 


Target 


Energy (GeV) 


Q 2 (GeV 2 ) 


X 


Completed experiments 


SLAC - E80, E130 


1976-1983 


e~ 


H-butanol 


^23 


1-10 


0.1-0.6 


SLAC - E142/3 


1992-1993 


e~ 


NH 3 , ND 3 


£ 30 


1-10 


0.03-0.8 


SLAC - E154/5 


1995-1999 


e~ 


NH 3 , 6 LiD, 3 He 


^ 50 


1-35 


0.01-0.8 


CERN - EMC 


1985 


M+ 


NH 3 


100, 190 


1-30 


0.01-0.5 


CERN - SMC 


1992-1996 


M + 


H/D-butanol, NH 3 


100, 190 


1-60 


0.004-0.5 


FNAL E581/E704 


1988-1997 


V 


P 


200 


~ 1 


0.1 <x F < 0.8 


Analyzing and/or Running 


DESY - HERMES 


1995-2007 


e + , e 


H, D, 3 He 


~ 30 


1-15 


0.02-0.7 


CERN - COMPASS 


2002-2012 


M + 


NH 3 , 6 LiD 


160, 200 


1-70 


0.003-0.6 


JLab6 - Hall A 


1999-2012 


e~ 


3 He 


< 6 


1-2.5 


0.1-0.6 


JLab6 - Hall B 


1999-2012 


e~ 


NH 3 , ND 3 


<6 


1.-5 


0.05-0.6 


RHIC - BRAHMS 


2002-2006 


P 


p (beam) 


2x (31-100) 


~ 1-6 


-0.6 < x F < 0.6 


RHIC - PHENIX, STAR 


2002+ 


V 


p (beam) 


2x (31-250) 


~ 1-400 


~ 0.02-0.4 


Approved future experiments (in preparation) 


CERN - COMPASS-II 


2014+ 




unpolarized H2 


160 


~ 1-15 


~ 0.005-0.2 








NH 3 


190 




-0.2 < x F < 0.8 


JLabl2 - HallA/B/C 


2014+ 


e~ 


HD, NH 3 , ND 3 , 3 He 


<12 


~ 1-10 


~ 0.05-0.8 



A. SLAC experiments 

SLAC experiments pioneered polarized DIS measure- 
ments and set many standards in polarized beam and 
target technologies. Their spin program focused on 
high statistics measurements of the inclusive asymme- 
tries. The first measurements of the proto n spin structure 



were performed by the e xperiments E 80 (Alguar det al. 



19761 . Il978h and E130 (|Baum et all . 1198(1 Il983h . foL 

lowed by a ser i es of high pre cision expe r iment s E142 



(lAnthonv et all Il996l). E14 3 (lAbe et all. Il998l) E154 



(|Abe et a/.l . [l997l) and E155 (jAnthonv et a?J . ll999ll2000h 
a decade later. These experiments utilized polarized elec- 
trons which were produced by laser photoemission and 
subsequently accelerated. The longitudinal polarization 
of the beam was frequently inverted and the polariza- 
tion measured using M0ller scattering. A rapid cycling 
of the beam and/or target polarization reduces system- 
atic uncertainties in the measured spin asymmetries re- 
lated to the stability of the experimental setup. Polarized 
target materials involved solid-state butanol and ammo- 
nia (NH3) for the proton and D-butanol, ND3 as well 
as 6 LiD for the deuteron. The target material, doped 
with a paramagnetic substance or irradiated with elec- 
tron beams, was polarized using dynamic nuclear polar- 
ization, which requires temperatures of about 1 K and 
strong magnetic holding fields. Such targets contain a 
considerable amount of non-polarizable nucleons, which 



is parametrized by the so-called dilution factor. This 
factor depends on all kinematic variables relevant for the 
process under study and needs, in principle, to be de- 
termined for each type of measurement. Typical values 
for polarized solid state targets range between 0.1 and 
0.2 with the exception of 6 LiD (0.4-0.5) and represent 
an important factor in the extraction of physical observ- 
ables from the measured ones. Information on the neu- 
tron structure was obtained either from the combination 
of measurements with proton and deuteron targets or 
by using a polarized 3 He target which is dominated by 
the neutron since the two proton spins in 3 He are anti- 
aligned. Here, polarization was obtained from optical 
pumping and adiabatic spin exchange. The target polar- 
ization was measured using the NMR technique. Scat- 
tered electrons were detected with magnetic spectrome- 
ters optimized for high-momentum-resolution and good 
electron identification. 



B. CERN experiments 

1. The EMC and SMC experiments 

Following the early measurements at SLAC, the Euro- 
pean Muon Collaboration (EMC) experiment performed 
at CERN in 1985 the first polarized DIS measurements at 
x < 0.1 down to x — 0.01 after a series of measurements 



(. 



of unpolarized nucleon and nuclear structure functions. 
The experiment used the polarized CERN muon beam up 
to momenta of 200 GeV and a solid-state irradiated am- 
monia target. Their low- a; measurements, accessible due 
to the high energy of the muons, suggested the break- 
down of the naive parton picture that quarks provide 



essentially all of the spin of the nucleon (jAshman et al 
19881 . Il989h . 



This triggered more detailed and precise measurements 
by the Spin Muon Collaboration (SMC) in 1992-1996, 
and by COMPASS (since 2002). The beam line and the 
principal ideas of the CERN muon experiments are de- 
scribed in the COMPASS Section 

After 1987 the focus was on the region x < 0.1 and 
the flavor-singlet axial-charge (Ellis-Jaffe sum-rule) for 
the neutron. The latter must deviate from the naive pre- 
diction in a similar way as for the proton in order to 
preserve the fundamental isovector Bjorken sum-rule for 
g\—gi- (These sum- rules are discussed below.) The SMC 
experiment could extend the measured x-range down to 
x = 0.004 (for Q 2 > 1 GeV 2 ) and established the va- 
lidity of the Bjorken sum-rule with measurements using 
polarized proton (but anol and ammonia ) and d euteron 
(D-butanol) targets (|Adeva et all Il993l Il998bl ). The 
large acceptance of the SMC spectrometer in the for- 
ward direction allowed them to present the first deter- 
minatio n of individu a l qua r k distr ibutions for different 



flavors (|Adeva et all Il996l Il998al) from semi-inclusive 
DIS. A dedicated polarimeter confirmed the validity of 
the bea m polarization obtained from Monte Ca rlo sim- 
ulations (Ada ms et all 2000; lAdeva et all Il994h used in 
the EMC, SMC and COMPASS analyses. 



2. The COMPASS experiment 



The COMPASS spectrometer (jAbbon et all 120071) is 



installed at the muon beam line of the CERN SPS accel- 
erator. A polarized muon beam of energy 160-200 GeV 
and with a polarization of about 80% impinges on a solid- 
state polarized target consisting of two or three cells with 
proton or deuteron target material polarized in oppo- 
site directions. The usable beam intensity is typically 
2 x 10 7 /s during a 9.6 s long spill. The repetition rate 
varies and is typically about 1/40 s. The muon polariza- 
tion arises naturally from the weak decay of the parent 
pions produced by the primary proton beam of 400 GeV. 
The momentum of each beam muon is measured in the 
beam momentum station. Downstream of the target, the 
scattered muon and produced hadrons are detected in 
a two-stage magnetic spectrometer with the two dipole 
magnets. 

Charged particles are tracked in the beam regions by 
scintillating fiber stations and by silicon detectors. In 
the inner region close to the beam, gaseous detectors of 
the micromegas and gas-electron-multiplier (GEM) types 



with high rate capabilities are deployed. The backbone 
of tracking in the intermediate region is multiwire pro- 
portional chambers (MWPCs). Finally, the large area 
tracking away from the beam region is covered by drift 
chambers and drift tubes. 

The velocity of charged particles is measured in a ring- 
imaging Cherenkov detector (RICH), which can separate 
pions and kaons from 9 GeV up to 50 GeV. The inner 
quarter of the photon detector is made of multianode- 
photomultiplicr tubes, while the outer part relies on MW- 
PCs with a photosensitive Csl cathode. 

The energy of charged particles is measured in sam- 
pling hadron calorimeters, while neutral particles, in par- 
ticular high-energy photons, are detected in electromag- 
netic calorimeters. They comprise lead glass modules as 
well as scintillator/lead "shashlik" modules in the inner 
high-radiation region. 

Event recording is triggered by the scattered muon, 
which is "identified" by its ability to traverse thick 
hadron absorbers, located just upstream of the Muon 
Wall detectors. The event selection is based on various 
systems of scintillator hodoscopes and logic modules ap- 
plying selection criteria like target pointing and energy 
loss in the scattering. The patterns causing a trigger 
were optimized by Monte Carlo simulations. The spec- 
trometer has about 250k read-out channels, which can be 
recorded with a frequency of 20 kHz for an event size of 
the order of 40 kByte. 

The heart of the experiment is the polarized target sys- 
tem. While the muon beam comes naturally polarized 
due to the parity violation in the decay of the parent 
pions, polarizing protons and deuterons is very difficult. 
Gas targets can not be used with the muon beam due 
to the low beam intensity compared to electron beams. 
An advantage of muon beams is the high muon energy, 
which presently can not be reached by electron beams. 
The polarized target system comprises a 2.5 T solenoid 
magnet, a 0.6 T dipole magnet, a 3 He/ 4 He dilution refrig- 
erator, a 70 GHz microwave system and an NMR system 
to measure the target polarization. The target material 
is cooled down to about 60 mK in frozen spin mode. The 
nucleons/nuclei are polarized by dynamic nuclear polar- 
ization which only is applicable for particular materials. 
In COMPASS irradiated ammonia (NH3) and lithium-6 
deuteride ( 6 LiD) were selected as proton and deuteron 
targets, respectively. Lithium-6 is very close to a system 
of a free deuteron and a helium-4 core and has essen- 
tially the same magnetic moment as the deuteron. Thus 
6 LiD corresponds to two deuterons plus a helium nucleus. 
Typically, polarizations of 85% for protons and 50% for 
deuterons were reached. A key feature of COMPASS is 
that both target polarizations are present simultaneously 
in separate target cells along the beam, e.g. "— >, <— " for 
the two-cell configuration until 2004 and "— », — >" for 
the three-cell configuration from 2006 onward. In the 
former configuration the length of the cells was twice 
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60 cm while in the latter it was 30 cm, 60 cm, 30 cm, 
respectively. Thus in an asymmetry measurement most 
systematic uncertainties cancel. Using the dipole and 
solenoid magnet, the magnetic field can be rotated from 
e.g. pointing downstream to transverse to upstream. The 
spin follows the magnetic field adiabatically and thus the 
spin orientations can be changed within 30 min. Such 
a held rotation is performed typically once per day for 
the longitudinal polarization in order to cancel poten- 
tially remaining systematic effects. The held can also be 
kept transverse for measurements with transverse target 
polarization. Here the polarization is inverted by repo- 
larizing typically once per week. In the shutdown year 
2005 the superconducting target magnet was replaced 
by a new one, increasing the angular acceptance from 
±70 mrad to ±180 mrad. 

The experiment is taking data since 2002. The years 
2008-2009 were dedicated to the hadron spectroscopy 
program of COMPASS with pion, kaon and proton 
beams. In 2012 the pion polarizability will be measured 
using a negative pion beam and a thin nickel target. A 
pilot run for deeply virtual Compton scattering and hard 
exclusive meson production is also planned for 2012. 



C. The HERMES experiment at DESY 

The HERMES experiment employed an innovative 
technique for the polarized target, which is very different 
from all other polarized DIS experiments. Gas targets of 
pure nuclear-polarized atoms of hydrogen or deuterium 
were used, which permit essentially background-free mea- 
surements from highly polarized nucleons with little or 
no dilution of the signal from unpolarized nucleons in 
the target. This choice eliminates one of the main sys- 
tematic sources in polarized DIS, the uncertainty in the 
determination of the dilution factor. 

The HERMES gas targets were highly longitudinally 
(~ 85%) or transversely (75%) polarized with the ability 
to invert the direction of the spin of the nucleons within 
milliseconds. Due to the low densities, however, such 
targets are only practicable in the high currents of storage 
rings. HERMES was operating from 1995 to 2007 at the 
HERA lepton storage ring, which provided electron or 
positron beams of typically 40 mA and with an energy 
of 27.5 GeV. In order to enhance the target density, the 
novel technique of a storage cell was used. Here, the 
gas was fed into a T-shaped open-ended elliptical cell 
coaxial to the lepton beam. The gas atoms underwent 
several hundred wall bounces before escaping from the 
ends where they were differentially pumped away by a 
large system of turbo-pumps. This increased the density 
by a factor of about 100 compared to free gas jet targets. 

The polarized atoms were injected into the cell from an 
atomic beam source based on Stern-Gerlach polarization 
filtering and radio-frequency transitions between atomic 



substates in a magnetic held (jAirapetian et all l2005a h 
A small sample gas diffused from the middle of the cell 
into a Breit-Rabi polarimeter which measured the atomic 
polarization, or into a target gas analyzer which mea- 
sured the atomic and the molecular content of the sam- 
ple. A magnet surrounding the storage cell provided a 
holding held defining the polarization axis and prevented 
spin relaxation via spin exchange or wall collisions. The 
cell temperature was kept at about 100 K, the value for 
which atomic recombination and spin relaxation during 
wall collisions are minimal. 

Stored high energy electron beams may become spon- 
taneously transversely polarized via a small polariza- 
tion asymmetry in the emission of synchrotron radi- 
ation by the beam particles as they are deflected by 
the magnetic fields of the rin g (Sokolov-Ternov ef- 
fect) (|Sokolov and Ternovl . 119641) . The beam polariza- 
tion grows and approaches asymptotically an equilibrium 
value with a time constant depending on the characteris- 
tics of the ring, for HERA typically 1/2 hour. Polariza- 
tions as large as 60% were achieved. Spin rotators and 
polarimeter s were essentia l components of the HERA lep- 



ton beam (Barber et all 1994 . 1993 ; Beckmann et al. 



2002 ; Buon and SteffenT 1986 ). Spin rotators in front of 



and behind the experiment provided longitudinal polar- 
ization at the interaction point and at one of the two 
beam polarimeters. The two beam polarimeters were 
based on Compton back-scattering of circularly polarized 
laser light. They continuously monitored the transverse 
and longitudinal polarization of the lepton beam. 

The HERMES spectrometer was designed to detect the 
scattered lepton and produced hadrons within a wide an- 
gular acceptance and with good momentum resolution. 
Particular emphasis was given to the particle identifica- 
tion capabilities which allowed for pion, kaon and proton 
separation over almost the whole momentum range. The 
HERA beam lines passed through the non-instrumented 
horizontal mid-plane of the spectrometer. A horizontal 
iron plate shielded the beam lines from the 1.5 Tm dipole 
field of the spectrometer magnet, thus dividing the spec- 
trometer in two identical halves. The geometrical accep- 
tance of ±170 mrad horizontally and ±(40 — 140) mrad 
vertically resulted in detected scattering angles ranging 
from 40 to 220 mrad. Tracking was provided by several 
stages of drift chambers before and after the spectrome- 
ter magnet. The combination of signals from a lead-glass 
calorimeter, a preshower detector, a transition radiation 
detector and a ring-imaging Cherenkov provided lepton 
identification with very high efficiency and purity bet- 
ter than 99% as well as pion, kaon and proton separation 
over almost the whole momentum range of 2-15 GeV. All 



compon ents are described in detail in lAckerstaff et al. 
(|l999af) . 



A recoil detector was installed in the target region for 
the last 1.5 years of HERMES data taking with unpolar- 
ized hydrogen and deuterium targets in order to enhance 
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access to hard exclusive processes, in particular to deeply 
virtual Compton scattering. 



D. JLab experiments 

Experiments at Jefferson Lab utilized the highest po- 
larization electron beams (85%) with energy ranging 
from 0.8 GeV close to 6 GeV. The technologies of polar- 
izing beam and target follow those pioneered and further 
developed at SLAC. 

The beam was provided by the Continuous Electron 
Beam Accelerator Facility (CEBAF), which used polar- 
ized electron guns based on a "superlattice" of a thin 
gallium arsenide (GaAs) layer on top of GaAs-phosphide 
bulk matter illuminated by circularly polarized photons 
from high intensity lasers. Subsequently, the polarized 
electrons passed up to five times the two linear accel- 
erators based on superconducting radio frequency tech- 
nology and connected by two recirculation arcs. The 
spin direction of the electrons was manipulated using 
the crossed electric and magnetic fields of Wien filters, 
which allow for rapid spin rotation. Their direction was 
inverted every about 30 ms. Beam polarimetry was em- 
ployed at several stages of the acceleration process. CE- 
BAF delivered polarized beams simultaneously to the 
three experimental halls (Hall A, B and C) with the op- 
tion to independently dial the energy and intensity. Typ- 
ical beam intensities ranged from a few nA in Hall B to 
over 100 /j,A in the other two halls. 

Longitudinal polarized solid state ammonia (NH3) tar- 
gets for the proton and ND3 for the deuteron were em- 
ployed at Hall B. These targets are based on similar tech- 
niques as discussed before for the SLAC and CERN ex- 
periments for both polarization and polarimetry. Hall 
A used a polarized 3 He target. The target polarization 
was measured by both the NMR technique of adiabatic 
fast passage and a technique based on electron param- 
agnetic resonance ([Romans et aLl . ll998l ). Average target 
polarizations of about 55% were obtained. 

Hall A and C were both instrumented with small accep- 
tance but high resolution spectrometers that could cope 
with the highest beam intensities but measured at fixed 
scattering angles. These spectrometers are equipped 
for high resolution tracking, precise time- of-flight mea- 



surem ents and lepton/hadron separation ([Alcorn et al. 
2004h . 

Hall B was instrumented with the CEBAF La r ge Ac 



ments, and electromagnetic calorimeters for neutral par- 
ticle identification. For charged particles, CLAS covered 
polar angles between 8° and 142° in the laboratory frame 
and between 60% and 80% of the azimuthal angles. 



E. Hadronic scattering experiments 

While deep- inelastic lepton-nucleon scattering has long 
been a standard tool of the trade in the study of unpo- 
larized and polarized nucleon structure, much has been 
learned from polarized hadronic scattering as well. The 
first high-energy primary polarized proton beams were 
achieved at the Zero-Gradient Synchrotron at Argonne 
National Laboratory in 1973. Proton beams there were 
initially accelerated to 6 GeV with a polarization of 
about 60%, and shortly thereafter polarized beams up 
to 12 GeV were achieved. In the 1990s at Fermilab, sec- 
ondary beams of polarized protons or antiprotons from 
lambda or antilambda decays opened up new kinematic 
regions for polarized hadronic scattering, with polar- 
ized beams of up to 200 GeV (^/s = 19 GeV). Polar- 
ized hadronic scattering experiments at center-of-mass 
energies more than an order of magnitude higher were 
achieved with the inauguration of the Relativistic Heavy 
Ion Collider for polarized protons in 2001. 



1. The Relativistic Heavy Ion Collider 

RHIC is located at Brookhaven National Laboratory 
in New York. RHIC was built to collide heavy ions at 
center-of-mass energies of up to 200 GeV per colliding 
nucleon pair and polarized protons at center-of-mass en- 
ergies ranging from 50 to 500 GeV. Collision of asym- 
metric species, i.e. different species in the two beams, is 
also possible due to independent rings with independent 
steering magnets. The first polarized proton collisions 
were achieved at a center-of-mass energy of 200 GeV in 
December 2001. 

The RHIC storage ring is 3.83 km in circumference 
and is designed with six interaction points (IPs) at which 
beam collisions are possible. Up to 112 particle bunches 
per ring can be injected, in which case the time between 
bunch crossings at the IPs is 106 ns. Polarizations of 
up to 65% for 100 GeV proton beams and about 60% 
for 255 GeV beams have been achieved. The maximum 

at 



ceptance Spectrometer (CLAS) (jMecking et all 120031) . y/s = 200 GeV and 2 x 10 32 cur 2 s" 1 at y/s = 510 GeV. 



The CLAS design was based on a toroidal field, generated 
by six superconducting coils arranged around the beam 
line. The six coils naturally divided the detector into six 
independent spectrometers, each of them containing a set 
of drift chambers for tracking, a gas Cerenkov counter for 
electron/pion separation, an array of scintillator counters 
for particle identification using time of flight measure- 



Three experiments have studied polarized proton 
collisions at RHIC. Th ere are two o n going large 



experimen ts, STAR (lAckermann et all [2003) and 



PHENIX (jAdcox et aLl . l2003l ). each of which have more 
than 500 collaborators total working on both the heavy 
ion and po larized proton p r ogram s, and the smaller 
BRAHMS (jAdamczvk et all |2003() experiment, with 
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fewer than 100 collaborators, which took data through 
2006. In additional to the program of proton spin struc- 
ture measurements, the transverse single-spin asymmetry 
in elastic proton-proton scattering has also been mea- 
sured to cons train the hadronic spin -flip amplitude in 



this reaction ([Adamczvk et all l2012bt ) . 



2. RHIC as a polarized p + p collider 

RHIC is the first and only high-energy polarized 
proton-proton collider in the world. A number of techno- 
logical developments and advances over the past several 
decades have made it possible to create a high-current 
polarized proton source, maintain the beam polarization 
throughout acceleration and storage, and obtain accurate 
measurements of the degree of beam polarization at vari- 
ous stages from the source to full-energy beams in RHIC. 
For an overview o f RHI C as a polarized-proton collider 
see lAlekseev et all (|2003l ) . In the case of polarized-proton 
running at RHIC, a pulse of polarized H~ ions from 
the source is accelerated to 200 MeV in the linac, then 
stripped of its electrons as it is injected and captured 
as a single bunch of polarized protons in the Booster, 
which accelerates the protons to 1.5 GeV. The bunch of 
polarized protons is then transferred to the Alternating 
Gradient Synchrotron (AGS) and accelerated to 24 GeV 
before injection into RHIC. Each bunch is accelerated in 
the AGS and injected into RHIC independently, with the 
two RHIC rings being filled one bunch at a time. The 
direction of the spin vector is selected for each bunch 
separately. The nominal fill duration is eight hours, af- 
ter which the beams are dumped and fresh beams are 
injected into RHIC. The bunch-by-bunch spin patterns 
in consecutive fills are varied in order to reduce potential 
systematic effects. 

Polarized proton injection uses an optically-pumpe d 
polarized H~ ion source (OPPIS) (IZelenski et al\ .E002). 
H~ polarization at the source of 85% has been achieved. 
There are two principal types of depolarizing resonances 
that may be encountered during acceleration, imperfec- 
tion and intrinsic resonances. Imperfection resonances 
refer to those due to errors in magnet currents or align- 
ments; intrinsic resonances are due to the (desired) beam 
focusing fields. A depolarizing resonance is encountered 
whenever the spin precession frequency is equal to the fre- 
quency with which a depolarizing field is crossed, leading 

to additive depolari zation effects. 

1978). 



Siberian snakes (Derbenev et al. 



spin-rotating dipoles, so named because of the beam 
trajectory through the magnets and the fact that they 
were developed at Novosibirsk in Russia, are used to 
overcome both imperfection and intrinsic resonances in 
RHIC. There are two snakes installed in each RHIC ring 
at diametrically opposite points along the rings. The two 
snakes in each ring rotate the spin vector 180° about per- 



pendicular horizontal axes, without perturbing the stable 
spin direction and with only local distortion of the beam 
orbit. In this way, all additive depolarizing effects from 
resonances are avoided. 

For RHIC to provide full-energy polarized beams, the 
polarization must be measurable at various stages of ac- 
celeration in order to identify and address possible origins 
of depolarization at each step. Only RHIC polarimc- 
try will be discussed here. There are two types of po- 
larimeters installed in RHIC. The proton-carbon (pC) 
polarimctcr takes advantage of a known analyzing power, 
A P N ss 0.01, in the elastic scattering of polarized protons 
with carbon atoms (p^ + C — > p^ + C) in order t o measure 



the beam polarization (jNakagawa et al . l2008h . This an- 
alyzing power originates from interference between elec- 
tromagnetic and hadronic elastic scattering amplitudes, 
which arc finite due to the anomalous magnetic moment 
of the proton. Calibration of the pC polarimeter to 
within an absolute beam polarization of less than 5% 
can be pr ovided by a po l arized hydrogen-jet-target po- 
larimeter ( Zelenski et al. . 2005). The technique is sim- 
ilar to that of the pC polarimeter in that the left-right 
transverse single spin asymmetry of elastically scattered 
protons due to the CNI process is measured. However, 
the analyzing power for this process of elastic proton- 
proton scattering at a given energy does not need to 
be known. Instead, the hydrogen-jet-target polarization, 
greater than 90%, is known to be tter than 2% in absolute 
polarization ( Okada et aZ.I . [2006) . The beam polarization 
can be obtained by measuring first the asymmetry due to 
the polarized target, averaging over the spin states of the 
beam, and then similarly measuring the asymmetry due 
to the polarized beam, averaging over spin states of the 
target, which are varied in time. The same (unknown) 
analyzing power, An, applies to both cases, allowing the 
beam polarization to be calculated. Low rates for this 
process mean that measurements performed over a long 
time are necessary. Because of this, the pG polarimeter, 
which can make measurements in less than ten seconds 
and provide immediate information on the stability or 
decay of the beam polarization from a few data points 
taken over the several hours of a fill, is essential. Results 
obtained from the two different polarimeters can later be 
compared in order to calibrate the pC polarimeter at any 
beam energy. 

The naturally stable spin direction through accelera- 
tion and storage in RHIC is transverse to the proton's 
momentum, in the vertical direction. Spin rotator dipole 
magnets hav e been used to achieve both radial and longi- 
tudinal spin (jMacKav et «Ll . l2003T ). The spin vector is ro- 
tated away from vertical immediately before the collision 
point and then back to vertical immediately afterwards. 
The rotators are located outside the interaction regions of 
the PHENIX and STAR experiments, giving both exper- 
iments the ability to choose independently whether they 
want longitudinally or transversely polarized collisions. 
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The BRAHMS experiment, having no spin rotators avail- 
able, focused on transverse spin measurements. The local 
nature of the spin rotator magnets means that the STAR 
and PHENIX experiments must each have their own way 
of checking the direction of the spin vector at their re- 
spective interaction regions. 

Observed azimuthal transverse single spin asymmetries 



in th e production of forward neutrons ([Bazilevskv et al 



20031 ) and forward charged particles can be used to pro- 
vide local polarimetry. These asymmetries are exploited 
by PHENIX and STAR to measure the degree to which 
the beam polarization is vertically transverse, radially 
transverse, or longitudinal. Maxima and minima of par- 
ticle production are expected at right angles in ip to 
the spin vector. For longitudinal polarization, no az- 
imuthal asymmetry is expected; a non-zero longitudinal 
single spin asymmetry would violate parity and is for- 
bidden in particle production via strong and electromag- 
netic processes. More information on local polarimetry at 



PHENIX an d STA R can be found in lAdare et all (|2007j) 



and lKirvlukl (120051) . 



3. RHIC experiments 

a. The PHENIX detector PHENIX was designed as a 
large, multi-purpose experiment with fast data acqui- 
sition and high granularity over a limited acceptance. 
There are two central arms with an acceptance cover- 
ing a pseudorapidity range | -^7 1 < 0.35 and A(p = | 
each in azimuth. The central arms include drift and pad 
chambers for momentum and position measurements, a 
ring-imaging Cherenkov detector primarily for electron 
identification, small-acceptance time-of-flight and aero- 
gel counters for charged hadron identification, and elec- 
tromagnetic calorimetry. Electronics-level triggering in 
the central arms uses information from the calorime- 
try and ring-imaging Cherenkov detector. There are 
two muon spectrometers covering a pseudorapidity of 
1.2 < \t]\ < 2.4, consisting of tracking chambers and 
muon identifier panels. Resistive plate chambers were 
added in 2011 and 2012 to improve triggering on high- 
momentum muons for W boson measurements. Forward 
electromagnetic calorimetry covering 3.0 < |?y| < 3.9 was 
added in 2006 and 2007, and silicon vertex detectors for 
heavy flavor measurements over I77I < 2.4 were added in 
2011 and 2012. 

For luminosity measurements, identical zero-degree 
hadronic calorimeters (ZDCs) are located in the RHIC 
tunnel at ±18 m from the nominal IP for all RHIC experi- 
ments. PHENIX also uses quartz Cherenkov beam-beam 
counters positioned around the beam pipe at ±1.44 m 
from the nominal interaction point as a minimum-bias 
trigger detector and for polarization-averaged as well as 
spin-dependent luminosity measurements. Collision rates 
for 500 GeV p + p running reach ^3 MHz, and the 



electronics-level triggers select events to reduce this rate 
to approximately 7 kHz of recorded data. 



b. The STAR detector The Solenoidal Tracker at RHIC 
(STAR) was designed as a large, multi-purpose detec- 
tor with wide acceptance, making it well suited for cor- 
relation measurements. The core of STAR is a time- 
projection chamber, which covers 2-7T in azimuth and 
has tracking capabilities over \rj\ < 1.3 and good parti- 
cle identification for \r)\ < 1. There is electromagnetic 
calorimetry for — 1 < 77 < 2. In the forward direc- 
tion, there is additional electromagnetic calorimetry for 
2.5 < rj < 4.0. Recent upgrades include a time-of-flight 
detector with 100 ps resolution for additional particle 
identification, and tracking based on GEM detectors for 
1 < r\ < 2 was partially installed for 2012 data-taking 
to enable charge-sign discrimination of forward electrons 
from W boson decays. 

In addition to the zero-degree hadronic calorimeters 
identical among the RHIC experiments, STAR has scin- 
tillator beam-beam counters positioned around the beam 
pipe covering 3.4 < |?7| < 5.0, which provide a minimum- 
bias trigger as well as spin-averaged and spin-dependent 
luminosity measurements along with the ZDCs. 



c. The BRAHMS detector The BRAHMS detector was a 
smaller experiment at RHIC designed for excellent mo- 
mentum measurement and charged particle identification 
over a very broad range of rapidities. It consisted of two 
movable spectrometer arms covering small solid angles, 
the Forward Spectrometer, which could be positioned as 
close as 2.3° from the beam pipe, and the Midrapidity 
Spectrometer, which could be moved to cover an angu- 
lar range from 30° < 6 < 95°. The spectrometer arms 
included five dipole magnets, time-projection chambers, 
multi-wire drift chambers, time-of-flight hodoscopes, and 
threshold as well as ring-imaging Cherenkov detectors. 
Global detectors consisted of a silicon array for multi- 
plicity measurements, threshold Cherenkov beam-beam 
counters for event vertex and timing determination as 
well as luminosity measurements, and ZDCs identical to 
those used by PHENIX and STAR. 



IV. THE PROTON SPIN PUZZLE 

We begin our discussion of physics results by first de- 
scribing how the small value of the "quark spin content" 
is obtained from polarized deep inelastic scattering 
and the first moment of the g\ spin structure function. 

In QCD the first moment of g\ is determined from the 
dispersion relation for polarized photon-nucleon scatter- 
ing and the light-cone operator product expansion. One 
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finds that the first moment of g\ is related to the scale- 
invariant axial charges of the target nucleon by 

/ dx g{(x,Q 2 ) 



12 9a ' 36 



^i 8) Vl + £ c NS^ s (Q)} 
/ e>i 

+^i 0) |inv{l + £c s ,af(g)}+O(^)+ 

£>1 ^ 

(12) 

Here g A \ g^ A and <7^|mv are the isovector, SU(3) 
octet and scale-invariant flavor-singlet axial-charges re- 
spectively. The flavor non-singlet cns£ and singlet cse 
Wilson coefficients are calculable in ^-loop perturba- 
tive QCD. These perturbative QC D coefficien t s hav e 
been calculated to O(o^) precision (jLarin et all . 1 19971 ) . 
For a s — 0.3 typical of the deep inelastic experiments 



one finds < 1 



v^3 i 



(Q)} 



= 0.85 and <^ 1 



S<Li c S£«f(Q)| = 0.96. The term ^ represents a pos- 
sible leading-twist subtraction constant from the circle at 
infinity when one closes th e contour in the complex plane 
in the dispersion relation ( Bassl . 2005 ). The subtraction 
constant affects just the first moment and corresponds to 
a contribution at Bjorken x equal to zero. 

In terms of the flavor-dependent axial-charges 

2Ms f _ l Aq = (p, s\qj^ 5 q\p, s) (13) 

the isovector, octet and singlet axial charges are 

g [ A ] = Ait - Ad 
g ( A ] = Au + Ad- 2As 
gf\i m /E{a s ) = gf = Au + Ad + As. (14) 

Here 



and g A are renormalization group invariants. We are 
free to choose the QCD coupling a s (p) at either a hard 
or a soft scale [i. The perturbative QCD expansion of 
E(a s ) remains close to one - even for large values of 
a s . If we take a s ~ 0.6 as typical of the infra-red then 
E{a s ) ~ 1-0.13-0.03 + .. . = 0.84+.. . where -0.13 and 
—0.03 are the 0(a s ) and 0(oq) corrections respectively. 

In the naive parton model g^ A is interpreted as the 
fraction of the proton's spin which is carried by the in- 
trinsic spin of its quark and anti-quark constituents. The 
experimental value of g^ A is obtained through measuring 
gi and combining the first moment integral in Eq.(12) 
with knowledge of g A and g A from other processes plus 
theoretical calculation of the perturbative QCD Wilson 
coefficients. 

The isovector axial-charge is measured independently 
in ne utron /3-decays (g A = 1.270+0.003 ( Beringer et all 
2012TH and the octet axial charge is commonly taken 
to be the value extracted from hyperon /3-decays as- 
suming a 2-parameter S U(3) fit (g^ = 0.58 ± 0.03 
(Clos e and Robertd . [199^ 1. However, it should be noted 



the uncertaint y quoted for g\ has been a matter o f 
some debate (Jaffc and Manoharl . Il990t iRatcliffd . 12004) . 
SU(3) symmetry may be badly broken and some have 
sugge sted that the error on effl should be as large as 
25% (jjaffe and Manoharl . Il99dh . A recent re-evaluation 
of the nucleon's axial-charges in the Cloudy Bag model 
taking into account the effect of the one-gluon-exchange 
hyperfine interaction and the pion cloud plus kaon loops 
led to the value g [ A ] = 0.46 ± 0.05 (|Bass and Thomasl . 
2010l ). The model reduction of g^ A from the SU(3) value 

(3) 

comes primarily from the pion cloud with g A taking its 
physical value. 

Deep inelastic measurements of g\ have been per- 
formed in experiments at CERN, DESY, JLab and 
SLAC. An overview of the world data on the nu- 
cleon's gi spin structure functio n is shown in Fig. 1 . 



E(a s ) — exp / da s 7(a s )//3(a s ) (15) This data is publi s hed in EMC ( Ashman et al. 



SMC (lAdeva et aLl . ll998bl) . E142 



is a renormalization group factor which corrects for the 
(two loop) non-zero anomalous dimension 7(oj s ) of the 
singlet axial-vector current 



Anthony et al. 



= uj^ 5 u + dj^ 5 d + S7 M 7 5 s 



(16) 



1989), 
1996), 



E143 (lAbe et all 119981). E 154 (lAbe et all Il997h . 
E155 ([Anthony et all |200d), E155 ( Anthony et al. 
19991). HERMES jAirapetian et all. l2007al). JLab 



which is close to one and which goes to one in the limit 



COMPASS ( Alekseev et all . l2010dt Alexakhin et al 



5 



( Dharmawardane et al . 2006 : IZheng et al. . 2004f). and 



Q 2 -> oo 
tion f3(a s 



The symbol /3 denotes the QCD beta func- 



11 - |/ )(a 2 /2n) + ... and 7 is given 



by j(a s ) = f {a s/ it) 2 + ... where / (=3) is the num- 
ber of active flavors ( Kodaira . 1980T) . The singlet ax- 
ial charge, g A |in V5 is independent of the renormaliza- 
tion scale fi and corresponds to g A \Q 2 ) evaluated in the 
limit Q 2 — > 00. The flavor non-singlet axial-charges g^ A 



20071 ). There is a general consistency among all data 
sets. The kinematic reach of the different experiments is 
visible in Fig. 2. COMPASS have the smallest- x data, 
down to x ~ 0.004. 

There are several striking features in the data. COM- 
PASS measurements of the deuteron spin structure func- 
tion gf show the remarkable feature that gf is consis- 
tent with zero in the s mall- z region between 0.004 and 
0.02 (jAlexakhin et all . 120071 ). In contrast, the isovec- 
tor part of g\ is observed to rise at small x as, g\ n ~ 



12 




FIG. 1 World data on xgi as a function of x for the proton 
(top), the deuteron (middle) and the neutron (bottom) at the 
Q 2 of the measurement. Only data points for Q 2 > 1 GeV 2 
and W > 2.5 GeV are shown. 



2 ,-o.22±o.o7 (jAlekseev et all l2010dl) and is much bigger 
than the isoscalar gf. This compares to the situation in 
the unpolarized structure function where the small- a; 
region is dominated by isoscalar gluonic exchanges. 



A. Spin sum-rules 

To test deep inelastic sum-rules it is necessary to have 
all data points at the same value of Q 2 . Next-to-leading 
order (NLO) QCD-motivated fits taking into account the 
scaling violations associated with perturbative QCD are 
used to evolve all the data points to the same Q 2 . First 
moment sum-rules are then evaluated by extrapolating 
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FIG. 2 World data for #i(x,Q 2 ) for the proton with Q 2 > 
1 GeV 2 and W > 2.5 GeV. For clarity a constant o, = 
0.28(11.6 — i) has been added to the gi values within a par- 
ticular x bin star ting; with i = for x = 0.006. (Also shown 
is the QCD fit of lLeader efaH lj2006l U 



these fits to x — and to i = 1, or using a Regge- 
motivated extrapolation of the data. Next-to-leading or- 
der (NLO) QCD-motivated fits discussed in Section V.C 
are used to extract from these scaling violations the par- 
ton distributions and in particular the gluon polarization. 

Polarized deep inelastic scattering experiments are in- 
terpreted in terms of a small value for the flavor-singlet 
axial-charge. For example, COMPASS found using the 
,(8) " 



SU(3) value for g\' (jAlexakhin el all 120071 ) and no lead- 
ing twist subtraction constant 



(0)| 

9A lpDIS,Q 2 



= 0.33 ± 0.03(stat.) ± 0.05(syst.). (17) 



(This deep inelastic quantity misses any contribution to 
^A^linv f rom a possible delta function at x = 0). When 
combined with g^ = 0.58 ± 0.03, the value of <7^| p dis 
in Eq.(17) corresponds to a negative strange-quark po- 
larization 



(°)i 

g \9 A IpDIS.Q 2 



9 f) 



-0.08 ± 0.01(stat.) ± 0.02(syst.) (18) 



- that is, polarized in the opposite direction to the spin 
of the proton. With this As, we obtain values for the up 
and down quark polarizations 

Auga^oo = 0.84 ± 0.01(stat.) ± 0.02(syst.) 
&d Q 2_> 00 = -0.43 ± 0.01(stat.) ± 0.02(syst.) (19) 
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FIG. 3 Convergence of the first moment integral of g\ as a 
function of the lower integration limit a; m i n for the Bjorken 
integral (isospin non-singlet) and the Ellis-Jaffe integral (iso- 
singlet) from the COMPASS proton and deuteron data at 
Q' 2 = 3 GeV 2 . The arrows indicate the theoretical expecta- 
tions. 



J 1 dxg\ +n are shown in Fig. 3. The Bjorken sum-rule 
and isosinglet integral converges in the measured x re- 
gion. Note that a large contribution, about 50%, of the 
Bjorken sum-rule comes from x < 0.15. The integral for 
the first moment of gi +n saturates at x ~ 0.05: the isos- 
inglet part of g± is close to zero in this measured range 
of small Bjorken x. 

The nucleon's second spin structure function gi is 
believed to satisfy the Burkhardt-Cottingham sum-rule 

dxg2 = ( Burkhardt and Cottingham . 1970h . The 
most precise measurements to date in polarized deep in- 
elastic scattering come from the SLAC El 55 and E143 
experiments, which report J ° 8 2 dx g v 2 = —0.042 ± 0.008 

for the proton and / ° 8 2 dx = —0.006 ± 0.011 for the 
deuteron at Q 2 = 5 GeV 2 (jAnthonv et al. 



20031) . E155 



estimate a contribution about 0.02 to the first moment 
of the proton 172 come from the x range between and 
0.02 from the twist-two (Wandzura-Wilczek) part of g^: 

s 2 w (*) = £f9i(y)-9i(z)- 



The non-zero value of Asg2_ i . 00 in Eq.( 18) is known as 
the v iolation of the Ellis-Jaffe sum-rule (jEllis and Jaffe , 
19741) . 

The extracted value of g^ | p diS required to be un- 
derstood by theory, and the corresponding polarized 

(8) 

strangeness, depend on the value of g\ '. If we in- 
stead use the value g^ = 0.46 ± 0.05 the correspond- 
ing experimental value of <7^| p dis would increase to 
5a 0) | p dis = 0.36 ± 0.03 ± 0.05 with 



As 



-0.03 ±0.03. 



(20) 



We shall discuss the value of As in more detail in Sections 
V and VI in connection with more direct measurements 
from semi-inclusive deep inelastic scattering plus global 
fits to spin data, models and recent lattice calculations 
with disconnected diagrams (quark sea contributions) in- 
cluded. 

The Bjorken sum- rule (|Biorkenl . Il966l Il97dh for the 
isovector part of g\ follows from current algebra and is 
a fundamental prediction of QCD. It has been confirmed 
in polarized deep inelastic scattering at the level of 5% 



£>1 



(<2)}- 



(21) 



B. Proton spin puzzles 

The results from polarized deep inelastic scattering 
pose the following questions: 

• How is the spin | of the proton built up from the 
spin and orbital angular momentum of the quarks 
and gluons inside? 

• Why is the quark spin content g^ | p dis so small? 

• How about g^' 7^ g^ ? What separates the values 
of the octet and singlet axial-charges? How reliable 
is the SU(3) value of gf 1 

• Is the proton spin puzzle a valence quark or 
sea/glue effect? 

• Can we extract information about the quark and 
gluon orbital angular momentum contributions 
from experiments, and with minimal model depen- 
dence? 

We next discuss the theoretical development and exper- 
imental work that has been performed to address these 
questions and the physics interpretation of present mea- 
surements. 



(3) 

The value of g A extracted from the most recent 
COMPASS data is 1 .28 ± 0.07(stat.) ± 0.010(syst.) 
([Alekseev et ali . l2010dl ) and compares well with the Par- 
ticle Data Group value 1.270 ±0 . 003 d educed from neu- 
tron beta-decays (jBeringer et aLl . l2012h . 

The evolution of the Bjorken integral dxg\~ n 
as a function of x m i n as well as the isosinglet integral big really is the OZI violation As 



C. Spin and the singlet axial-charge 

There are two key issues involved in understanding the 
small value of g^lpDis 1 the physics interpretation of the 
flavor-singlet axial-charge g^ and possible SU(3) break- 



ing in the extraction of g\ from hyperon /3-decays. How 

3 (9 A IpDIS 



,(°)| 
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Theoretical QCD analysis leads to the formula 



9X 



(0) 



(22) 



partons 



see 



Altarclli and Ross (1988), Ef remov and Tervaev 
1988[hlCarlitz et all (|l988h . iBass et all (|l99lh and lBass 
20051 ). Here Ag palt0 ns is the amount of spin carried 



by polarized gluon partons in t he polarized proton with 
a s Aq ^ constant as Q 2 — > o o ( Altarelli and Ross , 19881 : 
1988h ): Agp ar ton S measures the 



Efremov and Tervaev 



spin carried by quarks and anti-quarks carrying "soft" 

0(P 2 ,m 2 ) where P 2 is a 



transverse momentum 



u 2 



typical gluon virtuality in the nucleon and m is the 
light quark mass. The polarized gluon term is associ- 
ated with events in polarized deep inelastic scattering 
where the hard photon strikes a quark or anti-quark gen- 
erated from photon -gluon fusion and carrying k\ ~ Q 2 
( Carlitz et all 1988 ). It is associated with the QCD ax- 
ial anomaly in perturbative QCD. Coo denotes a poten- 
tial non-perturbative gluon top ological con tribution with 
support only at Bjorken x = (jBassl 2005). This term is 
discussed in Section VI on theoretical understanding. It 
is associated with the possible subtraction constant in the 
dispersion relation for g\ . If non-zero it would mean that 
lim e _j.o / dxgi will measure the difference of the singlet 
axial-charge and the subtraction constant contribution; 
that is, polarized deep inelastic scattering measures the 
combination g^ 0) | P Dis = 9^ - C oo- 

Possible explanations for the small value of <7^| p dis 
extracted from polarized deep inelastic experiments that 
have been suggested in the theoretical literature in- 
clude screening from positive gluon polarization, possible 
SU(3) breaking in the isosinglet axial-charges g^ and 
gjP, negative strangeness polarization in the nucleon, a 
possible topological contribution at x = plus c onnec - 
tions to axial U(l) dynamics discus sed in Fritzsch dl9891. 
Narison et al\ (Il995l) . IShorel (120081) and lBasd (jl999bl) . 

The two-loop QCD evolution factor E{a s ) in Eq.(15) 
is associated with the polarized gluon term which car- 
ries all the scale dependence. The quark spin contribu- 
tion Aq par tons and the subtraction constant in Eq.(22) are 



QCD scale invariant. The quark spin ter m Ag 



also known a s the J ET and chiral sche me (IChen 
Leader et al. 



1998) and AB scheme ( Ball et 




version of quark polarization - see Section VI. C. The 
growth in the gluon polarization Ag ~ l/ot s a t large Q 2 
is compensated by growth with opposite sign in the gluon 
orbital angular momentum. 

One would like to understand the dynamics which yield 
a small value of the singlet axial-charge extracted from 
polarized deep inelastic scattering and also the sum-rule 
for the longitudinal spin structure of the nucleon 



1 1 V- A 



Ag 



(23) 



where L q and L g denote the orbital angular momentum 
contributions. Operator definitions of the different terms 
or combinations of terms i n thi s equat i on are discussed in 
Jaffe and Manoharl (Il990l), [p dl997bh. IShore and White 



(2000), Bakke r et al. (12004) 



Bassl (I2005I), IChen et al 



2008), IWakamatsvJ feoiol) and lLeaderl (|201lh . We dis 



cuss orbital angular momentum and attempts to measure 
it in Sections VI- VII. 

There is presently a vigorous global program to dis- 
entangle the different contributions. Key experiments 
include semi-inclusive polarized deep inelastic scattering 
(COMPASS and HERMES) and polarized proton-proton 
collisions (PHENIX and STAR), as well as deeply vir- 
tual Compton scattering and hard exclusive meson pro- 
duction to learn about total angular momentum (COM- 
PASS, HERMES and JLab). Single spin observables in 
semi-inclusive scattering from transversely polarized tar- 
gets is sensitive to orbital angular momentum in the pro- 
ton. 



V. QUARK AND GLUON POLARIZATION FROM DATA 

Key observables needed to understand the small value 
of the singlet axial-charge g^lpDis are the polarized 
strangeness and polarized glue in the nucleon. The 
search for polarized strangeness has inspired a dedi- 
cated experimental program with semi-inclusive deep in- 
elastic scattering. Furt her, much activity was moti- 
vated by the discovery of Altarelli and Ross] ( 1988 ) and 
Efremov and Tervaev! ( 1988 ) that polarized glue makes a 
scaling contribution to the first moment of g±, a s Ag ~ 
constant. If there were a large negative contribution 
— 3^Ag with e.g. gluon polarization of the order of 
Ag ~ 2.5 at Q 2 = 10 GeV 2 , then this could reconcile 
the small measured value of IpDiS with the naive par- 
ton model expectation of about 0.6 through Eq. (22). 
This suggestion sparked a vigorous and ambitious pro- 
gram to measure Ag. Interesting channels include gluon 
mediated processes in semi-inclusive polarized deep in- 
elastic scattering (COMPASS and HERMES) and hard 
QCD processes in high-energy polarized proton-proton 
collisions at RHIC. 



A. Valence and sea polarization 

Semi-inclusive measurements of fast pions and kaons in 
the current fragmentation region with final state particle 
identification can be used to reconstruct the individual 
up, down and strange quark contributions to the pro- 
ton's spin. In contrast to inclusive polarized deep inelas- 
tic scattering where the g\ structure function is deduced 
by detecting only the scattered lepton, the detected par- 
ticles in the semi-inclusive experiments are high-energy 
(greater than 20% of the energy of the incident photon) 
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FIG. 4 Semi -inclusive longit udinal double-spin asymm etries for identified pions and kaons from COMPASS (jAlekseev et all 
l2009bll2010cT ) and HERMES l[ Airapetian et all l2OT5aT ) for the proton (left) and the deuteron (right) as function of x at the Q 2 
of the measurements. The error bars and bands indicate the statistical and systematic uncertainties, respectively. 



charged pions and kaons in coincidence with the scat- 
tered lepton. For large energy fraction z = Eh/E 1 1 
the most probable occurrence is that the detected w 



and K ± contain the struck quark or anti-quark in their 
valence Fock state. They therefore ac t as a tag of the 
flavor of the struck quark ( Close , 1979t ) . 



In leading order the virtual-photon-proton double-spin 
(cross-section) asymmetry is 



A lp (x, 



J2^ h elA q (x,Q 2 )jl ia dzD^z,Q 2 ) 
E q , h e 2 q q(x,Q 2 )Jl ia dzDh(z,Q^ 



(24) 



where z min ~ 0.2. Here D h q {z,Q 2 ) = J dp 2 (z , p 2 , Q 2 ) 
is the fragmentation function for the struck quark or anti- 
quark to produce a hadron h (= tt , K ± ) carrying energy 
fraction z = Ey l jE 1 in the target rest frame; Aq(x, Q 2 ) 
is the quark (or anti-quark) polarized parton distribution 
and e q is the quark charge. Note the integrati o n ove r 



1998al) and the HERM ES dAckerstaff etail Il999ri 



Airapetian et all 12004 l2005ai) experiments. The 



mo st recent measurements f rom HERMES are reported 
Airapetian et all (l2008ch and from COMPASS in 



1991) 



Alekseev et all (|2010c 

The experimental strategy has been to measure the 
asymmetries A\ for charged hadron production and sep- 
arated charged pion and kaon production from proton 
and deuteron targets. There is good agreement between 
the COMPASS and HERMES data in the kinematic re- 
gion of overlap - see Fig. 4. Flavor-separated polarized 
quark distributions for valence and sea quarks are then 
extracted from the data using fragmentation functions 
that have been fitted to previous hadron production data, 
with the most accurate taken to be those from the DSS 



group (|de Florian et all 120071 ) from a global fit to single- 



the transverse momentum p t (jClose and Milner 
Since pions and kaons have spin zero, the fragmentation 
functions are the same for both polarized and unpolar- 
ized leptoproduction. The fragmentation functions for 
u — » 7r + and d — > tt~ are known as "favored" (where the 
fragmenting quark has the same flavor as a valence quark 
in the final state hadron) ; the fragmentation functions for 
u — > 7r~ and d — > 7r + are known as "unfavored" . 

This program for polarized dee p inelastic s c atter- 
was pioneered by the SMC (|Adeva et all , Il996l l2005af) and SMC (|Adeva et al 



hadron production in e + e~, ep and pp collisions. 

The polarizations of the up and down quarks are pos- 
itive and negative respectively, while the extracted sea 
polarization data are consistent with zero - see Ta- 
ble II which includes measurements from COMPASS 



mg 



(Alek seev et all [2010c), 



HERMES ([Airapetian et al 
' 1998al) . 



J 



The COMPASS and HERMES determinations of the 
sum of strange and anti-strange polarisation As (a;) are 



r 



shown together in Fig. 5, plotted in the combination 
xAs(x). There is no evidence in the semi-inclusive data 



1G 



TABLE II First moments for valence quark and light-sea polarization from SMC, HERMES, and COMPASS. For each ex- 
periment the integrated sea is evaluated from data up to x — 0.3 and, for SMC, assuming an isospin symmetric polarized 
sea. 

Experiment x-range Q 2 (GeV 2 ) A^, Ack Am Ad 

SMC 0.003-0.7 10 0.73 ±0.10 ±0.07 -0.47 ±0.14 ±0.08 0.01 ±0.04 ± 0.03 0.01 ± 0.04 ± 0.03 

HERMES 0.023-0.6 2.5 0.60 ± 0.07 ± 0.04 -0.17 ± 0.07 ± 0.05 0.00 ± 0.04 ± 0.02 -0.05 ± 0.03 ± 0.01 
COMPASS 0.006-0.7 10 0.67 ± 0.03 ± 0.03 -0.28 ± 0.06 ± 0.03 0.02 ± 0.02 ± 0.01 -0.05 ± 0.03 ± 0.02 



in 

< 
x 



0.1 



COMPASS 
HERMES 




/o°oo4 dx (&u ~ Arf ) = 0.06 ± 0.04(stat.) ± 0.02(syst.). 
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FIG. 5 COMPASS (lAlekseev et all l2010ch and HERMES 
(jAirapetian et all l2008ct ) results for the strangeness polar- 
ization iAs(i) as function of x. The data are obtained in a 
leading-order analysis of SIDIS asymmetries (including those 
for charged kaons ) and u sing the DSS fragmentation functions 
(|de Florian et aUl2007D . 



for large negative strange quark polarization in the nu- 
cleon. The HERMES data covers the region 0.02 < 
x < 0.6, where the extracted As is consistent with a 
zero or small positive value. This data integrates to 
0.037 ± 0.019 ± 0.027 (|Airapetian et all . 



Cm dxAs 



2008c) in contrast with the negative value for polarized 
strangeness, Eq.(18), extracted from inc lusive measure- 



ments of q-i . COMPASS measurements (jAlekseev et al 



2009U l2010cf ) show no evidence of strangeness polar- 
ization in the region x > 0.004 with the integrated 
As = -0.02 ±0.02 ±0.02. 

The precise value of As extracted from semi-inclusive 
scattering may be affected by any possible future im- 
provement in the accuracy of the kaon fragmentation 
functions D^(z). However a drastic change in the ratio 

/ dzDf + 1 f dzD^ + would be needed to bring the first 
moment of As extracted from semi-inclusive scattering in 
agreement with the inclusive value, Eq.(18), obtained us- 
ing the SU(3) value of (jAlekseev et aZ.l . l2010cf) . More 



experimental data, especially on kaon fragmentation pro- 
cesses, are needed for improved precision on strangeness 
polarization in the nucleon. 

Semi-inclusive data are consistent with a small pos- 
itive or zero isospin asymmetry in the polarized sea 
Au - Ad. For the COMPASS data at 3 GeV 2 one finds 



For HERMES data at 2.5 GeV 2 f° ? 3 dx(Au - Ad) 



0.048 ± 0.057(stat.) ± 0.028(syst.) (|Airapetian et al 



2005a). These values compare with the unpolarized 

-0.118 ± 0.012 from 
A 



sea measurement f Q dx(u — d) 
the E866 experiment at FNAL (jTowell et all 1200 



compi lation of theoretical predictions is given in iPend 
(|2003l ). Meson cloud models predict small negative 



isosp i n asymmetries in the po l arized sea (|Cao and Signal 



2001 



, Kumano and Mivamal . 120021 ) whereas statistical 
(iBourrelv et all l2002h and chiral quark soliton models 
(|Wakamatsul . I2003T ) predict positive values. The COM- 



PASS and HERMES results are consistent with these pre- 
dictions within uncertainties. 

The VF-boson 
RHIC 



(Bunce et al. 



pro duction program at 
2000) will provide additional 



flavor-separated measurements of polarized up and down 
quarks and anti-quarks. At RHIC the polarization of 
the u, u, d, and d quarks in the proton is being measured 
directly using W boson production in ud — > W + and 
du — > W~ . The charged weak boson is produced through 
a pure V— A coupling and the chirality of the quark and 
anti-quark in the reaction is fixed. The W is observed 
through its leptonic decay W — > Iv, and the charged 
lepton is measured. Measurement of the flavor-separated 
anti-quark helicity distributions via W production in 
p ± p collisions is complementary to measurements via 
SIDIS in that there is no dependence on details of the 
fragmentation process, and the process scale, Q 2 « M 2 ^ 
is significantly higher than any data from existing fixed- 
target polarized DIS experiments. A parity-violating 
asymmetry for W + production in p + p collisions at 
y/s = 500 GeV consistent with predictions based on anti- 
quark helicity distributions extracted from SIDIS has 



alread y been obser ved by both PHENIX (jAdare et 



l2011bl ) and STAR (|Aggarwal et a/.l . l201l[ ) based on data 
collected in 2009. Considerably improved results are 
expected from data taken in 2011 and 2012 with higher 
luminosities and polarization. 

An independent measurement of the strange-quark 
axial-charge could be made through neutrino-proton elas- 
tic scattering. This process measures the combination 



(Au - Ad - As) 



0.01s. 



(0) 



where the small last 



term is a correcti on from heavy-quarks which has been 
calculated to LO (jKaplan and Manoharl 1198811 and NLO 



17 




FIG. 6 Jet production from quark-gluon scattering in 
polarized proton-proton collisions. 



( Bass et al. . 20021) accuracy. The axial-charge measured 
in vp elastic scattering is independent of any assumptions 
about possible SU(3) breaking, the presence or absence 
of a subtraction at infinity in the dispersion relation for 
<7i and the x ~ behav ior of q\ . 



In a recent analysis (|Pate et all 120081 ) of parity vio- 



lating quasi-clastic electron and neutrino scattering data 
between 0.45 and 1 GeV 2 (from the JLab experiments GO 
and HAPPEx and the Brookhaven experiment E734) , the 
axial form-factor was extrapolated to Q 2 = and favored 
negative or zero values of As with large uncertainty. 



B. Gluon polarization 

Polarized proton-proton scattering is sensitive to the 
ratio of polarized to unpolarized glue, Ag/g, via leading- 
order interactions of gluons, as illustrated in Fig. [6] The 
first experimental attempt to look at gluon polariza- 
tion was made by the FNAL E581/704 Collaboration us- 
ing a 200 GeV polarized proton beam and a polarized 
proton target. They measured a longitudinal double- 
spin asymmetry All for inclusive multi-7 and 7r°7r° pro- 
duction consistent with zero within their sensitivities, 
suggesting that Aq/q is not so larg e in the region of 



0.05 < x g < 0.35 (|Adams et aLl . 1199 



COMPASS was conceived to measure Ag via the study 
of the photon-gluon fusion process, as shown in Fig. 7. 
The cross-section for this process is directly related to 
the (polarized) gluon distribution at the Born level. 
The experimental techni que consists of the re construc- 



Pt hadrons ([Ageev et al 



tion of char med mesons ( Alekseev et all |2009cT > or high 



cess Ag. The high-py particle method leads to samples 
with larger statistics, but these have higher background 
contributions from QCD Compton processes and frag- 
mentation. High-pT hadron production was also used 
in ear ly attempts to ac c ess glu on polariza tion by HER- 



MES ([Airapetian et al. 
2004) and the most recent 



2000al) and SMC (|Adeva et al 
HERMES determination 



jAdolph et all l2012d) 



(Airapcti an et al\ | 2010cf) and COMPASS measurement 



These measurements in lepton-nucleon scattering are 
listed in Table [TTT1 for the ratio of the polarized to unpolar- 
ized glue Ag/g and shown in Fig.|S]for leading-order (LO) 
analyses of the data. The data cluster around x g ~ 0.1. 
There is no evidence in the data for non-zero gluon po- 
larization at this value of x a . 




FIG. 7 Production of a cc pair in polarized photon gluon 
fusion is being used to measure gluon polarization in the po- 
larized proton. 

The chance to measure Ag was a main physics drive 
for polarized RHIC. Experiments using the PHENIX 
and STAR detectors are investigating polarized glue in 
the proton. Measurements of Ag/g from RHIC are 
sensitive to gluon polarization in the range 0.02 < 
x g < 0.3 (Vs = 200 GeV) and 0.06 < x g < 0.4 
(y/s = 62 A GeV) for the neutr al pion All measured by 
PHENIX (jAdare et all l2009alibh and inclusive jet pro- 
duction measured by STAR a t 200 GeV center-of-mas s 



energy ([Abelev et all l2008bl : lAdamczvk et all l2012al) 



20061) in the final state to ac 



Additional channels s ensitive to Ag at RHIC have 
been publish ed as well ( Abelev et al. . 20091 : Adare et al. 
2011ail2012h . 



Preliminary results from PHENIX and STAR us- 
more recen t 200 G eV d ata than tho s e publ ished 



m g, more recen t ZOO i rev d ata tnan tno s e publ isned 
in Udare eiail (|2009bh and lAbelev et all (|2008bl ) are 
shown in Figs. [9] and [TO] In Fig. [9] the longitudinal 
double-spin asymmetry in neutral pion production mea- 
sured by PHENIX based on combined data from 2005, 



2006 , and 2009 is shown as a function of pr (jManion . 
201l[ ). The data are compared to a calculation using 
helicity distributions extracted from a global fit to po- 
larized world data from DIS, semi-in clusive PIS, and 
proton-prot on collisions ("DSSV") (|de Florian et al. 



20081 120091) and found to be consistent within the ex- 
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TABLE III Polarized gluon measurements from deep inelastic experiments. 

Experiment process (x g ) (p 2 ) (GeV 2 ) Ag/g 

HERMES (Airapetian et al. 2000a) hadron pairs 0.17 ~2 0.41 ± 0.18 ± 0.03 

HERMES CAirapetian et al. 2010c) inclusive hadrons 0.22 1.35 0.049 ± 0.034 ± 0.010±o.Jgg 

SMC fAdeva et al. 20041 hadron pairs 0.07 -0.20 ± 0.28 ± 0.10 

COMPASS fAgeev et al, 2006 ; Procureur. 2006) hadron pairs, Q 2 < 1 0.085 3 0.016 ± 0.058 ± 0.054 

COMPASS (Adolph et al. 2012c) hadron pairs, Q 2 > 1 0.09 3 0.125 ± 0.060 ± 0.063 

COMPASS (Alekseev et al. 2009c ; Franco. 2010) open charm (LP) 0.11 13 -0.08 ± 0.21 ± 0.11 
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FIG. 8 Gluon polarization Ag/g from leading-order analy- 
ses of hadron or hadron-pair production as function of the 
probed gl uon momentum fr action x g . Also shown a r e NLO 
fits from Ide Florian et al\ (I2009D and iLeader et al\ (|2010l ). 
(Figure taken from lAdolph etal\ i|2012d ).) 



perimental uncertainties; these data are also consistent 
with zero. Figure [TU] shows the longitudinal double-spin 
asymmetry in single-inclusive jet pr oduction measure d by 
STAR based on data taken in 2009 (|Diawothol . l201ll ) and 
provides the first evidence for non-zero gluon polarization 
in the proton. The results are compared to calculations 
based on the DSSV helicity distribution function extrac- 
tion as well as two earlier extractions based only on inclu- 
sive DIS data, GRSV-std and GRSV -zero correspondin g 



to Ag ~ 0.4 and 0.1 at Q 2 ~ 1 GeV 2 (jGlilck et al.l.l2001l) . 
See the following section for more details about fits to he- 
licity distributions. A given pt bin for single inclusive jet 
or hadron production generally samples a wide range of 
x g values. However, dijet measurements in p + p colli- 
sions provide better constraints on the x g values probed. 
Preliminary S TAR results for dijet production have also 
been released ( Walker , 20 111) and confirm the non-zero 
double-spin asymmetry seen in single jet production. 

While there is now evidence in the RHIC data that 
gluon polarization in the proton is non-zero, the mea- 
surements indicate that polarized glue, by itself, is not 
sufficient to resolve the difference between the small value 
of g^lpDis and the naive constituent quark model pre- 
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FIG. 9 The longitudinal double-spin asymmetry in 7r° pro- 
duction as measured by PHENIX in 200 GeV polarized pp 
collisions as a function of the transverse tv° momentum pr 
compared to a calculation u sing the DSSV polarized parton 
distributions l|Manionl . [20Tll ). 



diction, ~ 0.6, through the polarized glue term — 3^Ag. 
Note however that gluon polarization Ag ~ 0.2 — 0.3 
would still make a significant contribution to the spin of 
the proton in Eq.(23). 
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FIG. 10 The longitudinal double-spin asymmetry in inclusive 
jet production as measured by STAR in 200 GeV polarized 
pp collisions as a function of the transverse jet momentum pt 
compared to a calculations using the DSSV , GRSV-s t d, an d 
GRSV-zero polarized parton distributions ((Diawotho. i201lh . 



C. NLO QCD motivated fits to spin data 

Global NLO perturbative QCD analysis are performed 
on polarization data sets including both lepton-nucleon 
and proton-proton collision data. The aim is to ex- 
tract the polarized quark and gl uon parton dist ribu- 
tions. These analysis , starting from lBall et al. ( 1996T ) and 



Altarelli et all (| 19971 ) . frequently use DGLAP evolution 
and are performed in a given factorization scheme. This 
QCD fit approach has more recently been extended to a 
global analysis of data from polarized DIS, semi- inclusive 
polarized DIS and high-energy polariz ed proton-proton 



collisions (|de Florian et all [2008, 2009) 



The separation of g\ into "hard" and "soft" contribu- 
tions is not unique and depends on the choice of "factor- 
ization scheme" . For example, one might use a kinematic 
cut-off on the partons' transverse momentum squared 
(fef > fi 2 ) to define the factorization scheme and thus 
separate the hard and soft parts of the phase space for the 
photon-parton collision. The cut-off /i 2 is called the fac- 
torization scale. The coefficients in Eq.(8) have the per- 
turbative expansion C q = 6(1 — x) + §^f 9 (x, Q 2 /^ 2 ) and 
C 9 = §^f 9 (x, Q 2 /^ 2 ) where the strongest singularities in 
the functions f q and f 9 as x — > 1 are l n(l — x)/(l — x 



and l n(T — x) respectively, see e.g. lLampe and Reva 
(I2000TI . The deep inelastic structure functions are depen- 
dent on Q 2 and independent of the factorization scale /x 2 
and the "scheme" used to separate the 7* -parton cross- 
section into "hard" and "soft" contributions. 

Examples of different "schemes" used in the lit- 
erature are the modified minimal subtraction (MS) 



(|Bodwin and Qiul . Il990i: I't Hooft and Veltmanl . Il972h to 
regulate the mass singularities which arise in scatter - 
ing from massless partons, the "AB" (jBall et all 119961 ) 



( Leader et aUll998D schemes. 



19961 ) or "JET" 
In the MS scheme the po- 
larized gluon distribution does not contribute explicitly 
to the first moment of g\ . In the AB, CI and JET schemes 
on the other hand the polarized gluon (axial anomaly 
contribution) a s Ag docs contribute explicitly to the first 
moment since L dx = — ~ see the spin decom- 
position in Eq.(22). 

The /i 2 dependence of th e parton distributions is gi ven 
by the DGLAP equations (jAltarelli and Parisil Il977t) 



at 2ir 



f ^AP OT (*)A£(y,t) 

j x y y 

+ 2/ f 1 ^AP qg (-)Ag(y,t) 
Jx V V 

^ f ^AP gq ( X -)AX(y,t) 

Jx y v 

+ f d ^AP gg ( X -)Ag(y,t) 



(25) 

where AE(i,i) = ^ q Aq(x,t), t = ln/1 2 and / is 
the number of active flavors. The splitting func- 
tions Pij in Eq.( 25) have been calcu lated at lead- 
ing order by lAltarelli and Parisil (Il977l) and at next- 
to-leading order by IZiilstra and van Neerven ( 1994 ). 



Mertig and van Neerven (1996) and IVogelsand (jl996h . 

The largest uncertainties in the QCD fits are associ- 
ated with the ansatz chosen for the shape of the spin- 
dependent quark and gluon distributions at a given in- 
put scale. Further, the SU(3) value of gf ] (= 0.58 ± 
0.03) is assumed in present fits though no significant 
change in the x 2 quality of the fits should be ex- 
pected if one instead took a value of g^ with pos- 
sible 20% SU(3) breaking includedH The values for 
the quark and gluon spin contents (AE and Ag) ob- 
tained in recent NLO fits are listed in Ta b le IV with re- 
sults quote d from iBliimlein and Bottcherl (|2010l) (BB10: 
DIS data). iNocera et all (I2012D (NFRR12: DIS data), 
Leader et all (|2010|) (LLS10: DIS and SIDIS data) and 
Hirai and Kumand (120091) fAAC0 8: DIS and RHIC data) 
and lde Florian et all (I200rI . l2009h (DSSV08: DIS, SIDIS 
and proton-proton collision data). 

The most complete fits in terms of maximum included 
data are from the DSSV group, which take the SU(3) 

(8) 

value for g A ■ One finds need for a large negative con- 
tribution to As from small x, outside the measured x 
range when SIDIS data is included. The values obtained 
in this approach for J dxAs(x) are As — —0.057 
with a; m i n = and about —0.001 with x nlln = 0.001. 



1 We thank S. Taneja and R. Windmolders for discussion on this 
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TABLE IV First moments of th e polarized singlet-quark and gluon distributions at the scale 4 GeV 2 in the MS scheme; values 
quoted from lNocera et al\ (|2012t> . 





DSSV08 


BB10 


LSS10 


AAC08 


NFRR12 


AE(Q 2 ) 


0.25 ±0.02 


0.19 ±0.08 


0.21 ±0.03 


0.24 ±0.07 


0.31 ±0.10 


A 3 (Q 2 ) 


-0.10 ±0.16 


0.46 ± 0.43 


0.32 ±0.19 


0.63 ±0.19 


-0.2 ± 1.4 



That is, to reproduce the SU(3) value of the octet axial- 
charge, the negative polarized strangeness obtained from 
inclusive g\ measurements gets pushed into the unmea- 
sured small- a; range, x < 0.004. It is interesting here to 
note that, historically (before COMPASS, HERMES and 
RHIC Spin), the proton spin puzzle was assumed to be 
associated with strangeness/sea/glue polarization in the 
newly opened kinematics of EMC, SLAC and SMC, x be- 
tween 0.1 and 0.01. We now have accurate SIDIS mea- 



surements down to x ~ 0.004 which show no evidence 
for large sea/glue polarization effects. With the SIDIS 
measurements of As, either one needs SU(3) breaking in 
the octet axial-charge or strangeness/glue effects at very 
small x. Wit hout including the most rece nt data from 
2009 or later, Ide Florian etcdl (|2008t l2009h find a best- 
fit full first moment J* dxAg(x) = -0.084 at Q 2 = 10 
GeV 2 . With a very large Ax 2 /x 2 = 2% allowed range, 
the truncated first moment J dxAg = 0.013+o™4 
was obtained. With these errors Ag is still not precise. 



A recent attempt to extract polarized parton distribu- 
tions from inclusive polarized deep inelas tic data using 



neura l network techniques is reported in (jNocera et al 
20121 ). In this approach no assumption is made about 



the functional form of the input distributions, greatly 
reducing the primary (and notoriously difficult to quan- 
tify) systematic uncertainty on parton distribution fits. 
The neural network method has already been used 
quite successfully in the parametrization of the unpo- 
larized parton distribution s , with results published at 
NLO 



in 2010 (IBall et al 
20121) 



2010) and now NNLO 



2012 (BalLsLoU [2012). However, in the case of the ap- 
plication of the neural network method in the extraction 
of polarized parton distributions, thus far only inclusive 



polarized PIS data has been incor 
BB10 fit (Blumlein and Bottcher. 



porated, similar to the 
20101) . 



VI. THEORETICAL UNDERSTANDING 

In relativistic quark models some of the proton's spin 
is carried by quark orbital angular momentum. One has 
to take into account the four-component Dirac spinor for 



the quarks ip = -7= 



/ 

ia-fg 



where / and g are functions 



of the spatial coordinates and N is a normalization fac- 
tor. The lower component of the Dirac spinor is p-wave 
with intrinsic spin primarily pointing in the opposite di- 
rection to spin of the nucleon. In the MIT Bag model, 
where quarks are confined in an infinite square well po- 
tential with radius R, one finds the depolarization factor 
N 2 drr 2 (f 2 — ^g 2 ) = 0.65 for A q in the proton with 
all quarks in the (Is) ground-state (jJaffe and Manohar , 



1990h . That is, 35% of the proton's spin content is shifted 
into orbital angular momentum through the confinement 
potential. 

More detailed calculations of non-singlet axial-charges 
in relativistic constituent quark models are sensitive to 
the confinement potential, effective color-hyperfine inter- 
action, pion and kaon clouds plus additional wavefunc- 

tion corrections (associated with center of mass motion) 

(3) 

chosen to reproduce the physical value of g A . 

This physics was recen tly investigated within 
the Cloudy Bag model ( Bass and Thomasl . l2010t 



Mvhrer and Thomasi l2008t iThomasi 120081 ). The Cloudy 



Bag was designed to model confinement and sponta- 
neous chiral symmetry breaking, taking into account 
pion physics and the manifest breakdown of chiral 
symmetry at the bag surface in the MIT bag. If we 
wish to describe proton spin data including matrix 
elements of J^ 5 , J^ 5 and J m b, then we would like to 
know that the model versions of these currents satisfy 
the relevant Ward identities (the divergence equations 
for these currents). For the scale-invariant non-singlet 

axial-charges g A and g^ , corresponding to the matrix 
elements of partially conserved currents, the model is 
well designed to make a solid prediction. 

The effective color-hyperfine interaction has the quan- 
tum numbers of one-gluon exchange (OGE). In models of 
hadron spectroscopy this interaction plays an important 
role in the nucleon-A and S — A m ass difference s, as well 
as the nucleon magnetic moments ( CloseL Il979l ) and the 
spin a nd flavor dependence of p arton distribution func- 
tions ( Close and Thoniasl 1988 ). It shifts total angular- 
momentum between spin and orbital contributions and, 
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therefore, also contributes to model calculations of the 



octet axial-charges (jMvhrer and Thomasl . 119881 ). With 
OGE included (together with a phenomenological wave- 
function renormalization to ensure g A takes the physi- 
cal value) , the model is in very good agreement with the 
SU(3) fit to the nucleon and hyperon axial-charges ex- 
tracted from /3-decays with g A predicted to be around 
0.6. 

Next, the pion cloud induces SU(3) breaking in 
the nucleon's axial-charges. The pion cloud fur- 
ther shif ts intrinsic spin into orbital angular mo- 



mentum ( Schreiber and Tho mas. 1988; Tsush ima et 



19881 ). Including pion and kaon cloud corrections gives 



(8) 

the model result g A = 0.46 ± 0.05 (with the corre- 
sponding semi-classical singlet axial-charge or spin frac- 
tion being 0.42 ± . 07 bef ore inclusion of gluonic effects) 
(|Bass and Thomasl |2010|) . With this Cloudy Bag value 
for g A the corresponding experimental value of g^ A | p dis 



9 A IpDIS 



0.36 ± 0.03 ± 0.05, consid- 

pDIS — 



would increase to _ 

erably reducing the apparent OZI violation | (g A ^ 

(8) 

g A ) that one needs to explain. 

A recent lat tice calculation with disconnected dia- 
grams mcluded (|Bali et alikoi'A gave As = -0.02±0.01 
in the MS scheme at 7.4 GeV 2 . This value compares with 
the Cloudy Bag prediction As ~ —0.01 before gluonic 
degrees of freedom are included. These numbers are in 
good agreement with the values extracted from polarized 
SIDIS data by COMPASS and HERMES with the DSS 
fragmentation functions. 

Gluon polarization has been investigated in bag and 
light-cone models and in studies of heavy-quark axial- 
charges. The nucleon's charm-quark axial-charge was 
interpreted to give an estimate of gluon pola r izatio n 



l A f(™c)l < 0-3 with a s (m 2 c ) = 0.4 (|Bass et all 1201 If ). 



This upper bound corresponds to |3fjAg| < 0.06. Val- 
ues of Ag ~ 0.3 a nd ~ 0.5 at 1 GeV 2 were obtained in the 
MIT B ag model dChen and Jil . l2008h and in a light-cone 
model (jBrodskv et all 19951 ) respectively. These theo- 
retical values are consistent with the extractions of gluon 
polarization from COMPASS, HERMES and RHIC Spin 
data. 

To understand CooH deep inelastic sum- rules are de- 
rived using the operator product expansion and the dis- 
persion relation for deeply-virtual photon-nucleon scat- 
tering. Two important issues with the dispersion are 
the convergence of the first moment integral at the high- 
est energies and any contribution from closing the circle 
in the complex momentum plane. The subtraction con- 
stant, if finite, corresponds to a constant real term in the 
forward Compton scattering amplitude. It affects just 
the first moment integral and thus behaves like a S(x) 



term with support only at x = 0. A subtraction constant 
yields a finite correction to the sum-rule obtained from 
integrating only over finite n on-zero values of Bjorken x. 
One can show (| Bassl . 120051 ) that non-local gluon topo- 
logical structure requires consideration of a possible 8{x) 
subtraction constant. Whether it has finite value or not is 
sensitive to the real ization of a xial U(l) symmetry b reak- 
ing by instantons ( Crewther . 1978 ; 't Hooft . 19861) and 
the importance of topological structure in the proton. 
The QCD vacuum is a Bloch superposition of states char- 
acterized by non-vanishing topological winding number 
and non-trivial chiral properties. When we put a valence 
quark into this vacuum it can act as a source which po- 
larizes the QCD vacuum with net result that the spin 
"dissolves" . Some fraction of the spin of the constituent 
quark is shifted from moving partons into the vacuum at 
x = 0. This spin contribution becomes associated with 
non-local gluon topology with support only at Bjorken 
x = 0. 

Valuable information about the spin puzzle also follows 
from looking at the x dependence of g\ . The small value 
of g A ^ or "missing spin" is associated with a "collapse" in 
the isosinglet part of g\ to something close to zero instead 
of a valence-like rise at x less than about 0.05 - see e.g. 
the <7i data in Fig. 1 and the convergence of dxg\ +n 

and f dxgi~ n in Fig. 3. This isosinglet part is the sum 
of SU(3)-flavor singlet and octet contributions. If there 
were a large positive polarized gluon contribution to the 
proton's spin, this would act to drive the small a: part o f 
the singlet part of g\ negative ( Bass and Thomasl . 1 19931 ) 



- that is, acting in the opposite direction to any valence- 
like rise at small x. However, gluon polarization mea- 
surements constrain this spin contribution to be small in 
measured kinematics meaning that the sum of valence 
and sea quark contributions is suppressed at small x. 
Neither the SU(3) flavor-singlet nor octet contribution 
breaks free in the measured small x region. Hence, the 
suppression of g\ +n at small x should be either an isos- 
inglet effect or a delicate cancellation between octet and 
singlet contributions over an order of magnitude in small 
Bjorken x. 

The g\~ n data are consistent with quark model and 
perturbative QCD counting rules predictions in the va- 
lence region x > 0.2 ( Bassl . 1999a). The size of g A 
forces us to accept a large contribution from small x (a 
non-perturbative constraint) and the rise in g\~ n is in 
excellent agreement with t he pr e diction g^~ n ~ x~ ' 22 
of hard Regge exchange ( Bassl . 2007a ) - in particu- 
lar a possible a\ hard-pomeron cut involving the hard- 
pomeron which seems to play an i mportant role inun- 



polarized deep inelastic scattering ([Cudell et al , 1999) 



In the notation of Eq.(12): p x = -|Coo jl + J2e>i c Sl a i (<2)} ■ 



and in t he proton-proton total cross-sec tion measured 
at LHC ( Donnachie and Landshofl . 201 lh . (Soft) Regge 
theory predicts that the singlet term should behave as 
~ N In x in the small x limit, with the coefficient N to be 
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determined from experim ent (jBass and Landshofll . ll994 ; 
Close and Roberts! . 1994 ). From the data, this normal- 
ization seems to be close to zero. 

Where are we in our understanding of the spin struc- 
ture of the proton and the small value of g^ | p dis ? Mea- 
surements of valence, gluon and sea polarization suggest 
that the polarized glue term — 3^Ag and the strange 
quark contribution As part0 ns m Eq.(22) are unable to re- 
solve the small value of <7^| p dis- Two explanations are 
suggested within the theoretical and experimental uncer- 
tainties depending upon the magnitude of SU(3) break- 
ing in the nucleon and hyperon axial-charges. One is a 
value of g^jp ~ 0.5 plus an axial U(l) topological effect at 
x = associated with a finite subtraction constant in the 
<7i dispersion relation. The second is a much larger pion 
cloud reduction of g^ to a value ~ 0.4. Combining the 
theoretical error on the pion cloud chiral corrections em- 
braces both possibilities. The proton spin puzzle seems 
to be telling us about the interplay of valence quarks 
with chiral dynamics and the complex vacuum structure 
of QCD. 

Orbital angular momentum (OAM) in rclativistic 
quark models (for example the MIT and Cloudy bag 
models) without explicit gluon degrees of freedom has 
the usual interpretation of rclativistic quantum mechan- 
ics. For QCD dynamics the definition of OAM is more 
subtle because of the gauge covariant derivative, meaning 
that quark orbital momentum is in principle sensitive to 
the gluon fields in the nucleon that the quarks interact 
with. 

Going beyond spin and helicity to consider also orbital 
and total angular momentum, several operator decom- 
positions have been proposed. Starting from the relation 
between a ngular momentum and the energy-momentum 
tensor, [jj| (|l997bl ) takes 



J„= / <Fx x x T„ 



d 6 x 



ip^—ip + ip^x x (—iD)tp 



(fx xx(ExB) 



(26) 



The gauge covariant derivative means that L q is a pri- 
ori sensitive to gluonic degrees of freedom. The J q and 
J g quantities here are amenable to QCD lattice calcula- 
tions and, in principle, measurable through deeply virtual 
Compton scattering. In an alternative approach, taking 
the + light cone component o f the QCD angular momen - 
tum tensor in A+ — gauge I Jaffe and Manohar 
proposed the operator decomposition 

q+\xxio\ q+ 



M+ 12 = -g|7sg+ + ffl I x x id 



where the gluon term in the gauge covariant derivative is 
no longer present through the gauge fixing. 

The connection between the quark and gluon total an- 
gular momentum contributions J q and J 9 and the QCD 
energy-momentum tensor all ows us to write down their 
LO QCD evolution equations ( Ji et al. , 19961) . The quark 



and gluon total angular momenta in the infinite scal- 
ing limit are given by J q (oo) = |{3//(16 + 3/)} and 
J g (oo) = ^{16/(16 + 3/)}, with / the number of active 
flavors - that is, the same scaling limit as the quark and 
gluon momentum contributions at infinite Q 2 . The Ji 
and Jaffe-Manohar definitions of orbital angular momen- 
tum satisfy the same (LO) QCD evolution equation, so, 
at LO, are equal in a model calculation if the glue con- 
tribution can be set equal to zero at a low-energy input 
scale. 

To obtain information about the quark "orbital angu- 
lar momentum" L q we need to subtract the value of the 



"intrinsic spin" S q 



^Aq measured in polarized deep 



inelastic scattering from the total quark angular momen- 
tum J q . This means that L q is scheme dependent with 
different schemes corresponding to different physics con- 
tent depending on how the scheme handles information 
about the axial anomaly, large-fc t physics and any possi- 
ble "subtraction at infinity" in the dispersion relation for 
g\. The quark total angular momentum J q is anomaly 
free in QCD so that axial anomaly effects occur with 
equal magnitude and opposite sign in L q and S q . When 
looking at physical observables that are sensitive to OAM 
and quark spin (with possible axial-anomaly contribu- 
tion) it will be important to identify which OAM defini- 
tion and which scheme quantity is most relevant to the 
observable - for example, in SIDIS the largest k t events 
are included in the MS version of Aq whereas they are 
omitted in the JET scheme version (lBassLl2003h . 

To connect quark model predictions with lattice cal- 
culations and fits to data it is necessary to use QCD 
evolution of the model results from the low-energy scale 
where the model applies up to the hard scale of deep 
inelastic scattering. Model calculations (and also lat- 
tice calculations without disconnected diagrams) of Aq 
are commonly understood to refer to the scale invari- 
ant version of this quantity, e. g. the chiral/ JET or 
AB scheme quark spin contributions in Eq. (22). One 
chooses a model ansatz for the gluon polarization and 
total angular momentum, typically Ag = J g = at the 
model input scale. For illustration, Fig. 11 shows the 
evolution of total and orbital angular momentum con- 
tributions in the Cloudy Bag from the model scale up 
to Q 2 = 4 GeV 2 . Various phenomenological investi- 
gatio ns (jMattinglv and Stevenso 3, 11994 Stcffens et al 



-2TrF +j ix x idjA^ + Tre+- y F +i A J (27) 



Il995l ) found that by going beyond leading order QCD 
(and including pions in the nucleon wavefunction) , the 
optimal fit to high-energy scattering data involved taking 
the running coupling a s about 0.6-0.8 at the low energy 
input scale. For this range of a s the scale dependence of 
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FIG. 11 Calculation of the NLO QCD evolution of J u , L d , 
Jd, L- g in the Cloudy Bag with model input scale Qo = 0.4 
GeV (|Thomas et all I^Olof ) . 



AE (in full QCD) through Eq.(15) converges well. (Go- 
ing to higher orders in the model fits to data and putting 
in pions raises the model input scale /io needed for the cal- 
culatio ns.) Table V compar es the results of lattice calcu- 



lations (Hagler et al 



20081 ) for up- and down-quark spin 
and total angular momentum with the Cloudy Bag model 
results and the values extracted from QCD fits to hard 
exclu sive reaction data, GPDs ([Goloskokov and Kroll , 
2009), and transverse single spin asymmetries, TMDs 
( Bacchetta and Radici , 201 lh . The lattice calculation in- 
volves connected diagrams only (no axial-anomaly con- 
tribution) plus chiral extrapolation. The QCD fit num- 
bers are central values modulo (possibly large) system- 
atic errors from the model functional forms of distribu- 
tions used in the fits. There is good convergence of the 
different theoretical values with "data". Here, one has 
L u ~ — Ld ~ 15% at the scale of typical deep inelastic 
measurements . 



In an alternative approach, the proton spin puzzle has 
also been addressed in the Skyrme model, where baryons 
emerge as topological solitons in the meson fields at large 
number of colors N c , and in the chiral quark soliton 
model (ChQSM), where explicit quark degrees of free- 
dom are also present in the model. The nucleon's axial 
charges in these models are sensitive to which mesons are 
included in the model and the relative contribution of a 
quark source and pure meson component. In an early cal- 
culation iBrodskv et al. found that g^ vanishes in 



a particular version of the Skyrme model with just pseu- 
doscalar mesons. Finite values of As ~ —0.08, close to 
the value obtained from inclusive g± measurements with 
good SU(3) assum ed for g A are found in the ChQSM 
(|Wakamatsul . l2007t ). 



TABLE V Model, lattice and fit extractions of angular mo- 
mentum contributions in the proton, quoted for 4 GeV 2 (ex- 
cept GPD at 2 GeV 2 ). 





Cloudy Bag 


Lattice 


GPD 


TMD 


Au 


0.85 ±0.06 


0.82 ±0.07 






Ad 


-0.42 ±0.06 


-0.41 ±0.07 






J u 


0.30 


0.24 ±0.05 


0.24 


0.24 


Jd 


-0.04 


0.00 ± 0.05 


0.02 


0.02 



VII. TRANSVERSE NUCLEON STRUCTURE AND 
ORBITAL ANGULAR MOMENTUM 

Confinement induces transverse hadronic scales in the 
nucleon with accompanying finite quark orbital angular 
momentum and finite spin-orbit couplings which can be 
probed in experiments. The search for orbital angular 
momentum has motivated new theoretical and experi- 
mental investigations of the three dimensional structure 
of the nucleon. Key observables in deeply virtual Comp- 
ton scattering (DVCS) and transverse single spin asym- 
metries in lepton-nucleon and proton-proton scattering, 
and also the large x limit of the down quark helicity dis- 
tribution are sensitive to orbital angular momentum in 
the nucleon. 

Hard exclusive reactions such as DVCS are described 
theoretically using the formalism of generalized parton 
distributions (GPDs) and probe the three-dime nsiona l 
spatial struc t ure of the n ucleon as re viewed inljil ( 1998 1. 
Goeke et all (|200lh and iDiehll (|2003t ). Ji has derived a 
sum-rule connecting the forward limit of GPDs to infor- 
mation about the quar k and g luon total angular momen- 
tum in the proton (0, 1997bl ). Considerable experimen- 



tal and theoretical effort has and continues to be invested 
aimed at accessing this information. 

Studies of single spin asymmetries for semi-inclusive 
meson production in high-energy lepton-nucleon and 
proton-proton collisions are sensitive to possible spin- 
orbit coupling both in the nucleon and in the final- 
state hadronization process; for a recent review see 



Barone et all (|2010aj ). One studies correlations between 
the transverse momentum (orbital motion) of partons, 
their spin and the spin polarization of the nucleon. 
The theoretical tools are transverse momentum depen- 
dent (TMD) distributions and fragmentation functions. 
TMDs probe the three-dimensional transverse momen- 
tum structure of the nucleon and are associated, in part, 
with finite orbital angular momentum. 

Experimental studies of three-dimensional nucleon 
structure have been pioneered at HERMES and JLab for 
GPDs and at COMPASS, HERMES and RHIC for TMDs 
in single spin asymmetry measurements. There has also 
been considerable theoretical effort aimed at model and 
lattice calculations of these observables. 
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In the rest of this Section we present the theory and 
present status of these new GPDs and TMD distribu- 
tions plus spin-orbit coupling in fragmentation and the 
prospects for future experiments including key observ- 
ables that will be studied. The aim for experiments 
should be to focus on observables that have the clean- 
est theoretical interpretation with minimal model depen- 
dence. 

Quark orbital angular momentum in the nucleon may 
also be manifest in future measurements of the large x 
behavior of the polarized down-quark distribution Ad/d 
and in the ratio of the proton's spin-flip Pauli form- factor 
to the Dirac form- factor at large Q 2 . These observables 
can be studied with the 12 GeV upgrade of JLab. Va- 
lence Fock states with non-zero orbital angular momen- 
tum induce a logarithmic correction to the QCD count- 
ing rules predictions for these observables. Perturbative 
QCD calculations which take into account orbital angular 
momentum give 



F 2 /Fx~(log 2 Q 2 /A 2 )/Q 2 



(28) 



for the ratio of Pauli to Dirac form-factors at large Q 2 
(iBelitskv et all l2003h. Form-factor measurements at 
JLab (|Gavou et all 12001; I Jones et all . |2000|) are consis- 



tent with this behavior and also with F2/F1 ~ 1 / \[Q 2 
for Q 2 between 4 and 6 GeV 2 , in contrast to the count- 
ing rules predict ion without orbital a n gular momentum 



19801 ). One also 



F 2 /F\ ~ 1/Q 2 (|Lepage and Brodskv 
finds a logarithmic correction to the leading large- a; be- 
havior of the negative-helicity spi n-dependent quark dis - 



tributions — (1 - a;) 5 log (1 - x) (jAvakian et all 120071) 



An interesting prediction here is that Ad/d should cross 
zero and become positive at a value x ~ 0.75 when this 
term is included, in contrast to the model expectation 
that crossing occurs at x ~ 0.5 when this orbital an- 
gular momentum effect is neglected. An accurate mea- 
surement of Ad/d at x close to unity would be very in- 
teresting if this quantity can be extracted free of uncer- 
taintie s from nuclear effects (|Kulagin and Melnitchouk . 
2008alibh in the neutron structure functions measured 
from deuteron or 3 He targets. 



A. Generalized parton distributions 

Observables in deeply virtual Compton scattering and 
deeply virtual meson production are sensitive to infor- 
mation about total angular momentum in the nucleon. 
In these hard exclusive reactions a deeply virtual pho- 
ton impacts on a nucleon target and a real photon or 
a meson is liberated from the struck nucleon into the 
final state, leaving the target nucleon intact. These pro- 
cesses can be described using the formalism of general- 
ized parton distributions (GPDs), involving the Fourier 
transfo r ms of off-diagonal nucleon matrix eleme nts (|jl . 
1997al lbl: lMueller et a/.l . [l^^lRadvusnkmlll997l) . 



The important kinematic variables are the virtuality of 
the hard photon Q 2 , the momenta p— A/2 of the incident 
proton and p + A/2 of the outgoing proton, the invari- 
ant four- momentum transferred to the target t = A 2 , the 
average nucleon momentum P, the generalized Bjorken 
variable k + = xP + and the light-cone momentum trans- 
ferred to the target proton £ = — A + /2p + . In the Bjorken 
limit, £ is related to Bjorken xb via £ = ib/(2 - sg). 
The generalized parton distributions are defined as the 
light-cone Fourier transform of the point-split matrix el- 
ement[f| 



1 



cP+t 



y tf\Mv)Mo)\p)v+=v±=o 

H{x,t,t)u(p') 1 +u(p) 



-(757 )a/9 



+E(x,Z,t)u(p / )a+^u(p) 

H(x,£,i)S(p')7 + 7 5 u(p) 
A+ 

+E(x^,t)u(p') l5 —u(p) 



(29) 

The physical interpretation of the generalized parton 
distributions (before worrying about possible renormal- 
ization effects and higher order corrections) is the follow- 
ing. Expanding out the quark field operators in Eq. (|29|) 
in terms of light-cone quantized creation and annihilation 
operators one finds that for x > £ (x < £) the GPD is 
the amplitude to take a quark (anti-quark) of momen- 
tum k — A/2 out of the proton and reinsert a quark 
(anti-quark) of momentum fc + A/2 into the proton some 
distance along the light-cone to reform the recoiling pro- 
ton. In this region the GPD is a simple generalization 
of the usual parton distributions studied in inclusive and 
semi-inclusive scattering which are for mally defined via 
light-cone correlation functions - see e.g. iBassI (|2005ri . In 
the remaining region — £ < x < £ the GPD involves tak- 
ing out (or inserting) a qq pair with momentum k — A/2 
and — k — A/2 (or k + A/2 and — k + A/2) respectively. 
Note that the GPDs are interpreted as probability am- 
plitudes rather than densities. The non- forward matrix 
elements give access to transverse degrees of freedom in 
the nucleon. 

In the forward limit the GPDs H and H are related to 
the parton distributions studied in deep inelastic scatter- 
ing 



H{x,C,t)\ s=t =o = q(x) 
H(x,£,t)\t=t=o = &q{x) 



(30) 



3 We work in the light-cone gauge A+ = (so the path-ordered 
gauge-link needed for gauge-invariance in the correlation function 
becomes trivial and set equal to one). 
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whereas the GPDs E and E have no such analogue. Inte- 
grating over x the first moments of the GPDs are related 
to the nucleon form-factors 
-l 



dxH(x, £, t) 



1 

+1 



dxE{x, £, t) 



Fi(t) 



i 

+i 



dxH(x,£,t) = G A (t) 



i 

+i 



dxE(x,£,t) = G P (t). 



(31) 



Here F\ and F2 are the Dirac and Pauli form-factors of 
the nucleon, and G a and Gp are the axial and induced- 
pseudoscalar form-factors respectively. (The dependence 
on £ drops out after integration over x.) 

GPDs contain vital information about quark t otal an - 
gular momentum in the nucleon. Ji's sum- rule Jjl ll997bh 



relates J q to the forward limit of the second moment in 
x of the spin-independent quark GPDs 



dxx 



H q (x,£,t = 0)+E g (x,£,t = 0) 



(32) 



The gluon "total angular momentum" could then be ob- 
tained through the equation 



1 

j a = - 2 . 



(33) 



In principle, it could be extracted from precision 
measurements of the Q 2 dependence of DVCS at 
next-to-leading-order accuracy where the quark GPDs 
mix with glue under QCD evolution or via J g = 
i j+l dxx{H3(x, £_,t = 0)+ E3(x, £, t = 0)} if the gluon 
GPD can be accurately measured in more direct ex- 
periments. In these equations J q and J g are defined 
through the proton matrix elements of the angular mo- 
mentum operators in Eq.(26). If information about J q 
can be extracted from experiments, then the correspond- 
ing quark "orbital angular momentum" can be deduced 
by subtracting the value of the quark spin content Aq 
extracted from deep inelastic scattering and polarized 
proton-proton collisions. 

Experimental attempts to access J q via Eq.(32) require 
accurate determination of the two unpolarized GPDs H 
and E. Measurements from a proton target are more sen- 
sitive to J u , the total angular momentum carried by up 
quarks, while the neutron (via a deuteron or 3 He target) 
is most sensitive to Jd- The experiments require high 
luminosity to measure the small exclusive cross-section, 
plus measurements over a wide range of kinematics in 
Q 2 , x and t (since sum- rule tests and evaluations depend 
on making reliable extrapolations into unmeasured kine- 
matics). One also needs reliable theoretical technology 
to extract the GPDs from the measured cross-sections. 





DVCS Bethe-Heitler 

FIG. 12 The leading DVCS and Bethe-Heitler processes. 



GPDs appear in the amplitudes for DVCS and hard ex- 
clusive meson production as convolutions with the hard 
scattering coefficient and only these so-called Compton 
form factors (CFF) are experimentally accessible. Vary- 
ing the photon or meson in the final state gives access 
to different flavor combinations of GPDs, like in semi- 
inclusive DIS. However, meson production is more sensi- 
tive to QCD radiative corrections and power corrections 
in 1/Q, and reliable theoretical description requires larger 
values of Q 2 compared to DVCS. Channels particularly 
sensitive to gluons in the proton are hard exclusive vec- 
tor meson production where both quark and gluon GPDs 
appear at lowest order in the strong coupling constant. 
There is a challenging program to disentangle the GPDs 
from the formalism and to undo the convolution inte- 
grals which relate the GPDs to measured cross-sections. 
In practice, the approach used is to constrain models of 
GPDs against experimental data in measured kinemat- 
ics. These models are then integrated to obtain the Ji 
moments of J u and Jd, which may then be compared to 
the predictions of QCD inspired model plus lattice cal- 
culations - see Table V. 

In the rest of this discussion we focus on deeply virtual 
Compton scattering. 



1. Deeply virtual Compton scattering 

Measurements of hard exclusive processes are much 
more challenging than traditional inclusive and semi- 
inclusive scattering experiments. These exclusive pro- 
cesses require a difficult full reconstruction of final state 
particles and their cross-sections are usually small, de- 
manding high luminosity machines. 

DVCS experiments have to be careful to choose the 
kinematics so as not to be saturated by a large Bethe- 
Heitler (BH) background where the emitted real photon 
is radiated from the incident lepton rather than from the 
proton target - see Fig. 12. 

Most of the DVCS program so far has focused on the 
DVCS-BH interference term. Use of different combina- 
tions of beam and target polarization plus changing the 
electric charge of the incident lepton beam gives maxi- 
mum access to most combinations of DVCS observables. 
Measurement of the DVCS-BH interference term - see 
Eq.(34) below - allows one to measure not only the size 
of the DVCS amplitude but also its phase - that is, it 
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gives separate information about the real and imaginary 
parts of the Compton form factors. 

Pioneering measurements of DVCS have been per- 
formed at DESY (HERMES, HI and ZEUS) and JLab 
(Hall A and Hall B), which complement each other in the 
covered kinematic phase space and the extracted observ- 
ables. 

The experiments use different measurement techniques 
to access exclusive reactions. The HERA collider experi- 
ments at DESY, HI and ZEUS, as well as CLAS (Hall B 
at JLab) have the advantage of nearly hermetic spectrom- 
eters, whereas the fixed target experiments HERMES 
and JLab Hall A had to deal with the restrictions caused 
by incomplete event reconstruction due to their forward 
spectrometers. Hall A and HERMES successfully em- 
ployed the so-called missing mass t echnique together 



with careful background sub traction (jAirapetian et al. 



2008b; ICamacho et all 120061 ). For Hall A the low beam 



energy and high resolution spectrometer allowed one to 
resolve pure elastic scattering from associated produc- 
tion with an excited nucleon in the final state. The 
latter contribution was treated as part of the signal in 
HERMES results. Very recently, beam-spin asymme- 
tries for a pu re DVCS sample have a lso been reported 
by HERMES (lAirapetian et all 120123 ) . In the fixed tar- 
get experiments the spin-dependent DVCS cross-sections 
have been explored using longitudinally polarized lep- 
ton beams with longitudinally (JLab and HERMES) 
and transversely (HERMES) polarized targets. HER- 
MES also took advantage of the available different beam 
charges. 

JLab experiments focus on kinematics dominated by 
valence quarks. Data from Hall A suggest leading twist- 
2 dominance of DVCS even a t the relatively low Q 2 of 



1.5-2.3 GeV 2 ( Camacho et all 12006) 



The HERA collider experiments HI and ZEUS mea- 
sured the DVCS cross-section close to the forward direc- 
tion with £ < 10 -2 , integrated over its azimuthal depen- 
dence, in an xb range where two-gluon exchange plays a 
major role in addition to the leading order quark-photon 
scattering process. Figure [13] shows the cross-section 
differential in t for differ ent ranges in Q 2 measured 



by H I (jAaron et all 120081 ) and ZEUS (|Chekanov et al. 
2009). The data are well described by the exponential 
behavior da/dt oc e~ f, '*L The distribution of partons in 
the transverse plane is then obtained from this depen- 
dence by a Fourier transform with respect to At (the 
transverse momentum shift in Eq. (29)) F(b, x. Q 2 ) o c 
/ d 2 A T exp- lbAT V 'dajdt (iBurkardtl l20o3 iDiehll . XMvh . 
The impact parameter provides an estimate of the trans- 
verse extension of the partons probed during the hard 
process. While DVCS data provide information about 
the transverse distribution of quarks in the proton, data 
on exclusive heavy vector meson production ( J/ty or T) 
describe the transverse distribution of glue at specific val- 
ues of x. 




FIG. 13 The t dependence of the DVCS cross-section for sev- 
eral values of Q 2 as measured by HI and ZEUS. The curves 
are results of fits of the form e -6 '* with b being related to the 
transverse extension of partons in the proton at a given x and 
Q 2 (see text). 



Th e full DVCS cross-section reads ([Diehl and Sapeta 
20051) 



da(£p —¥ ijp) 



da uu 



da, 



DVCS 

uu 
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PiSt da LT 



e £ da\ju 
egPiS L da[ L 
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eiP e dai v 
egSh da\j L 

epSx dafjrp + St da^^ 



DVCS 
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P( St da^T < ~'^ 
Pzda^cs 

S L daB v L CS 



DVCS 



(34) 



Here the first subscript U, L on da indicates an unpo- 
larized or longitudinally polarized lepton beam and the 
second subscript U, L, T denotes an unpolarized, longi- 
tudinally or transversely polarized proton target; P( is 
the lepton beam polarization; Sl and St denote longi- 
tudinal and transverse proton polarization. The various 
cross-section terms depend on the azimuthal angle cj) be- 
tween the lepton scattering plane and the photon pro- 
duction plane, and, in case of a transversely polarized 
proton target, also on the azimuthal angle (j)s between 
the lepton plane and the transverse target spin vector. 
Equation ([54]) indicates the large variety of observables 
accessible with polarized beams and/or targets. 

Of particular interest is also the dependence on the sign 
of the charge of the beam lepton ei , which allows one to 
disentangle contributions from the pure interference term 



and the DVC S term as pioneered in lAirapetian et al. 
(l2007bLl2008bl) . 
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TABLE VI Linear combinations of Compton form factors 
(CFF) in the DVCS-BH interference terms. Here, Fi and 
F 2 are the electromagnetic form factors. Subleading terms 
not shown are suppressed in a wide range of kinematics. 



Target polarization 



CFF combination 



FiH + + F 2 )H - jKF 2 E 



unpolarized / charge 

longitudinal FiH + £(Fi + F 2 )H 

transverse oc sin(^ — (f>s) F 2 H — Fi E + . . . 

transverse oc cos(0 — (j>s) F 2 H — Fi ^E + . . . 



As an example for the azimuthal dependence of the 
cross-section we give the expression for the interference 
term fo r the case of an unpo larized target and polarized 



beam ( Belit skv et all l2002f) 



Joe 



' 3 



ci cos(n(f>) + A 4 sin(7i0) ] . (35) 



The proportionality involves a kinematic factor and the 
lepton propagators of the BH process; A is the helicity of 
the incoming lepton. The Fourier coefficients provide 
an experimental constraint on the real part of the Comp- 
ton form factor and on the imaginary part. Their 
relation to linear combinations of Compton form factors 
and hence to the respective GPDs is listed in Table IVT1 
A specific Fourier coefficient can be accessed experimen- 
tally by weighting the cross-section with the respective 
azimuthal modulation. 

The DVCS-BH interference term was extracted by 
varying the electric charge of the incident lepton (HER- 
MES) and studying polarization observables, varying the 
beam or target helicity (JLab and HERMES). JLab ex- 
periments have focused on studying their accessible ob- 
servables fully differentially. HERMES explored the ad- 
vantages of using simultaneously polarization and charge 
observables to cleanly isolate the interference term and 
obtained the most complete set of DVCS observables 
measured so far providing access to all interference terms 
listed in Eq. (j34f . 

Figure M shows a summary of the HERMES DVCS 
measurements with polarized proton a nd deuterium tar- 
gets at their avera ge kinem at ics (lAirapetian et al 



2008bl l2009cL l2010blldl . l2011allbl . l2012bj|d). Here, A c is 
the charge asymmetry and Axy ar e the polarization- 
dependent asymmetries with X and Y indicating the 
beam and target polarization, respectively, which could 
be longitudinal (L) or transverse (T). The subscript I in- 
dicates an extraction of the pure interference term. The 
measured asymmetries are subject to a harmonic expan- 
sion with respect to the azimuthal angle(s) as given by 
the superscript of Axy in the Figure. These data de- 
noted by squares in Fig. 14 show results extracted from 
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FIG. 14 Overview of all DVCS azimuthal asymmetry am- 
plitudes measured at HERMES with proton and deuterium 
targets, given at the average kinematics. 



a DVCS samp le with kinemat i cally c omplete event re- 
construction (jAirapetian et al\ . \2012a ). The dependence 
on the kinematic variables t, Q 2 , and Xb was explored 
for each observable. 

An example of the high statistics data from JLab is 
shown in Fig. [15] for the beam-spin asy mmetry Aj,\j 
meas ured fully differentially by CLAS (|Girod et all . 
20081 ). The presented data contain an admixture of the 



4 sm0 and /l 8111 ?' 
^LUJ tmu A LU,DVCS 



contributions from the interfer- 
ence and pure DVCS terms, which cannot be separated 
here. CLAS also provides measurements of A s ^\f and 
20061 ). 



A s ™ 2 ^ (jChen et al 



2. The quest for orbital angular momentum and GPD 
parametrizations 

Of the two GPDs H and E entering Ji's sum-rule, 
Eq.(32), measurements with unpolarized targets but 
longitudinally polarized beams and also beam-charge 
asymmetries are mainly sensitive to H. As indicated 
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FIG. 15 The leading beam-spin asymmetry amplitude 
a(t) = differential in t, x and Q 2 as measured 

by C LAS dGirod et all [2008). An earlier CLAS measure- 
ment (Stcpanva n et aZ.L I2001T) is indicated by the red square. 
The op en green triangles repre sent the cross-section data from 
Hall A [|Camacho et all I2006T ). 



in Table IVI1 transverse target polarization provides 
kinematics-wise unsuppressed access to E. The GPD 
E is essentially unknown. In contrast to H, it is not 
related to a deep inelastic parton distribution in the 
forward limit; E describes helicity flip at the proton 
vertex and requires finite orbital angular momentum in 
the nucleon. Contributions from E to most DVCS ob- 
servables are damped by kinematic factors ~ \t\/Mp, 
with the average \t\ value generally much smaller than 1 
GeV 2 in the experiments. To access E requires DVCS 
and/or vector meson production asymmetry measure- 
ments with transversely polarized nucleon targets. It 
may also be accessed through the beam polarization 
dependence of DVCS with a neutron target because 
of th e different size of th e form-factors for the neu- 
tron (jBelitsky et all 120021 ). Measurements have been 



performed al r eady for all c hanne ls (jAirapetian et al. 



Jd + Jtt/5.0 = 0.18 ± 0.14(expt.) (|Mazouz et all 120071 ) 
whereas HERMES results from the proton gave J„ + 



2008bl ) in 



J d /2.8 = 0.49 ± 0.17(expt.) (jAirapetian et 
the same model. 

To go further and perform global fits of GPDs to hard 
exclusive observables one faces several challenging theo- 
retical issues. Parameterizations of GPDs have to deal 
with two longitudinal variables instead of one plus the t 
dependence of DVCS. It is also not yet known whether 
relatively simple and smooth functions like those used 
in QCD fits to deep inelastic data are sufficient to de- 
scribe GPDs. A reliable parameterization of GPDs might 
therefore require a larger number of moments than em- 
ployed in usual QCD parton descriptions. In addition, 
the dependence of the functions on the variable x is 
not directly accessible, as x represents a mute variable 
which is integrated over. In the interpretation of DVCS 
observables one has to deal with complex amplitudes; 
the GPDs are embedded in the Compton form factors 
which relate to the measured cross-sections. Despite 
these complications and the early stage of global fitting 
for GPDs, many results have been obtained in recent 
years fitting to different hard exclusive scattering data. 
Int erested readers are re fe rred t o the original literature 
in Vanderhaeghen et al. (ll999T). Goloskokov and Kroll 



(2008), Guidal (2010) 



Kumericki and Muellerl (|2010l ) and 



Goldstein et ali (|2011l ). This phenomenology is comple- 



mented by progress in latt i ce QCD calculations o f GPD 



moments ( Bromm el et 
Hagler et all 120081) . 



2007; Gockcler et al.. 2007 



B. Transversity, transverse-momentum-dependent 
distributions and fragmentation functions 

Striking single spin asymmetries associated with spin- 
momentum correlations (expected with parton orbital 
angular momentum) were first observed in the 1970s. 
Using a 12 GeV polarized proton beam from the Ar- 
gonne National Laboratory Zero Gradient Synchrotron 
on a fixed target, up to 40% more positive pions 
were produced left of the beam when the beam was 
polarized up, and up to 20% more ne gative pions 



were produced to the right of the beam (jKlem et 
19761 ) . These measurements were confirmed by simi 



2008bl l2009at iMazouz et 'all 120071) . Despite the lack of lar experiments dAntille et aLl.ll98dlA~pokin et art ll99fi 



precision for these observables, attempts to extract in- 
formation about quark total angular momentum have 
been performed by fitting the oretical models o f GPDs 



Dragoset et al. . 1978t ISaroff et aZfl990 ). but it was not 



to the DVCS meas urements ( Airapetian et al . 2008ri 



Mazouz et al . 120071 ). Although this analysis is very 
model dependent, the results agree (surprisingly) well 
with model and lattice expectations, e.g. the calcu- 
lations reported in Table V. For e xample, within the 
model of lVanderhaeghen et all (|l999h JLab Hall A DVCS 
measurements from the neutron were interpreted to give 



until the 1990s that a theoretical framework was devel- 
oped to attempt interpreting them. 

Single-spin asymmetries have now also been observed 
in proton-proton collisions at RHIC, where they reach 
up to ~ 40%, and in lepton-nucleon collisions at COM- 
PASS, HERMES and JLab, where they are typically 
5 — 10%. Single-spin asymmetries for hadron produc- 
tion from transversely polarized targets tell us about 
spin-orbit coupling in the nucleon and/or in the fragmen- 
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tation process. Transverse-momentum-dependent distri- 
butions simultaneously describe the dependence on lon- 
gitudinal momentum fraction of the parton within the 
parent hadron as well as the parton's transverse momen- 
tum. Similarly, transverse-momentum-dependent frag- 
mentation functions describe the dependence on longitu- 
dinal momentum fraction of the produced hadron with 
respect to the scattering parton as well as the hadronic 
transverse momentum with respect to the jet axis. 

We next introduce these distributions and frag- 
mentation functions, and discuss their phenomenol- 
ogy. In QCD there are eight leading-twist quark 
T MDs. These are listed in Ta ble IVIII and discussed 
Mulders and Tangerman ( 1996t and Bacchetta et al. 



TABLE VII Leading-twist transverse-momentum dependent 
parton distributions. U, L, and T stand for unpolarized, 
longitudinally polarized, and transversely polarized nucleons 
(rows) and quarks (columns) respectively. 



The three distributions highlighted in boldface 
survive integration over transverse momentum k t . These 
yield the unpolarized parton distribution, fi(x,k t ), the 
spin-dependent parton distribution gi(x,k t ) and the 
trans versity distribution hi(x,k t ). The other five dis- 
tributions do not survive integration over k t . They de- 
scribe correlations between the quark transverse momen- 
tum with the spin of the quark and/or the spin of the 
parent nucleon, viz. spin-orbit correlations. The three 
TMDs denoted by h describe the distribution of trans- 
versely polarized partons. They are chiral-odd distribu- 
tions and appear only in observables involving two chiral- 
odd partners, such as Drell-Yan processes (two chiral-odd 
parton distributions) or SIDIS (chiral-odd parton distri- 
bution and the Collins fragmentation function discussed 
below). The three distributions (the Sivers distri- 
bution), (the Boer-Mulders distribution) and h^ T 
(pretzelosity) require orbital angular momentum in the 
nucleon since they involve a transition between initial 
and final nucleon states whose orbital angular momen- 
tum differ by ALf = ±1 (Sivers and Boer-Mulders) or 
ALf = ±2 (pretzelosity). The "worm-gear" functions 
h^ L and g^ T link two perpendicular spin directions and 
are also connected to quark orbital motion inside nucle- 
ons. 
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Transverse momentum distributions have been stud- 
ied most in semi-inclusive DIS experiments where they 
appear in combination with the usual unpolarized frag- 
mentation function D{z,p t ) or, in case of the chiral-odd 
TMD distributions, with a chiral-odd Collins fragmen- 
tation function H^-(z,p t ) discussed in Section VII. B. 2. 
One measures the azimuthal distribution of the produced 
final-state hadron with respect to the virtual-photon axis. 
Each species of TMD comes with a different angular mod- 
ulation in the semi-inclusive cross-section allowing it to 
be projected out to yield infor mation about the differ- 
ent spin-momentum correlations (jBacchetta et aLl . l2004l) . 
All these leading-twist TMDs have been measured in 
semi-inclusive DIS over the last decade. However, sev- 
eral have been the focus of more intense studies and we 
focus on those here. 

The different modulation combinations are listed in 
Table VIII together with present experimental measure- 
ments. Results quoted at = 18 GeV are from COM- 
PASS, 7.4 GeV from HERMES and 3.5 GeV from JLab. 
Here <fi is the angle between the lepton direction and the 
plane spanned by the exchanged photon and tagged final- 
state hadron, e.g. a high-energy meson; 4>s is the angle 
between the lepton direction and the transverse nucleon 
target spin. The convolution is taken over the involved 
transverse momenta of the quark and the hadron pro- 
duced in the fragmentation process. 



When one projects out the terms with different az- 
imuthal angular dependence summarized in Table VIII, 
the COMPASS, HERMES and JLab data suggest that 
the Sivers, Collins and Boer-Mulders effects are all 
present in the proton target data - see below. JLab data 
from CLAS reveal a clear signal for the worm-gear- 1 dis- 
tribution; there is some hint for a non-zero worm-gear-2 
distribution (with low significance) and (so far) no sig- 
nificant signal for prctzolosity in the proton. For the 
deuteron target, there is evidence for a Boer-Mulders ef- 
fect from COMPASS and HERMES. The Sivers, Collins, 
worm-gear and pretzelosity effects are all consistent with 
zero in the deuteron target data. The Collins and Sivers 



effects observed in the proton data therefore contain a 
predominant isovector contribution. 

We next focus on the Sivers, Boer-Mulders and Collins 
effects. 



1. The Sivers and Boer-Mulders TMD distributions 



The Sivers distribution was first proposed in Siversl 
(Il990h in an attempt to explain the large transverse sin- 
gle spin asymmetries observed in the 1970s and 1980s. It 
describes the correlation between the transverse momen- 
tum k t of the struck quark and the spin 5* and momentum 
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TABLE VIII Experimental access to the leading twist TMD distributions in SIDIS with unpolarized (U), longitudinally (L) 
or transversely polarized (T) beam (modulation first subscript) and/or target (modulation second subscript). 



Modulation Combination 




Target 


Observed 


Measurement 


Distribution name 


GeV 


type 


hadron types 




sin(0 + (f)s)uT hi <g> Hi 


18 
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h ± ,it ± ,K ± ,K° 


Ageev et al. (2007) and Alekseev et al. (2009a) 


Transversity 




P 




Adoloh et al. (2012a) and Alekseev et al. (2010b) 






P 


7V ± ,K ± 


Drelim. Pesaro (2011) 




7.4 
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AiraDetian et al. (2005b. 2010a) 
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n ± 


Oian et al. (2011) 
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18 
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h ± ,7T ± ,K ± ,K° 


Ageev et al. (2007) and Alekseev et al. (2009a) 
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P 
P 


7V ± ,K ± 


Adoloh et al. (2012b) and Alekseev et al. (2010b) 
Drelim. Pesaro (2011) 
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AiraDetian et al. (2005b. 2009b) 
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Oian et al. (2011) 
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prelim. Sbrizzai (2011) 
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Airaoctian a/. (2012a) 
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u 


7T + 


Osiocnko et al. (2009) 
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Drelim. Kotzinian (2007) 
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h ± 


prelim. Parsamyan (2011) 
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prelim. Paooalardo (2010) 
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Alekseev et al. (2010a) 
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Airapetian et al. (2003) 
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Avakian et al. (2010) 
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prelim. Kotzinian (2007) 
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prelim . Parsamyan (2011) 
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prelim. PaDDalardo and Diefenthaler (2011) 






3.5 
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Oian (2012) 



p of its parent nucleon 



f q / p f(x, k t ) = ff(x, k't) - fi t q {x, h t ) 



S-(ktxp) 



Al 



The k t dependence means that the Sivers distribution is 
sensitive to non-zero parton orbital angular momentum 
in the nucleon, though the mapping from Sivers observ- 
ables to quark (and gluon) orbital angular momentum is 
(so far) model dependent with present theoretical tech- 
nology. 

The Sivers distribution has the interesting property 
that it is odd under time reversal. Due to this fea- 
ture, such a correlation was believed to be forbidd en 
for more than a decade. Then Brodskv et al. ( 2002al fb) 
showed that, with initial- or final-state interactions, the 
Sivers effect could be non-zero in QCD processes. Final- 
state interactions in SIDIS can generate the azimuthal 
asymmetry before the qua r k fra gments into hadrons. 
Shortly afterwards, ICollinsI (|2002l) realized that initial- 
state color interactions in the case of Drell-Yan and 
final-state interactions in the case of SIDIS would lead 
to a process-dependent sign difference in the Sivers dis- 



tribution. SIDIS measurement s (lAdolph et al. 



(36) lPesard . l2011tlQian et al 



Airapetian et all [2005b , l2009ri lAlekseev et al. 



2012b; 



2010b; 



20111 ) , suggest sizable asymme- 



tries at the level of about 5 — 10% for a proton and a neu- 
tron target, while Drell-Yan measurements are planned 
for the future. 

A qualitative picture of the Sivers distribution can al- 
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FIG. 16 Sivers amplitudes for c harged pions mea s ured b y 
HERMES with a proton target (|Airapetian et all l2009bl ). 
The Sivers amplitudes for K + (not shown here) appear to 
be nearly twice as large as those for 7r + . 
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ready be deduced from SIDIS measurements. The non- 
zero amplitudes shown in Figs. [16] and [17] were ob- 
tained with a proton target. For HERMES the ampli- 
tude includes a kinematic factor depending on the ratio 
of transverse-to-longitudinal photon flux, which in the 
COMPASS data is divided out. Since scattering off u 
quarks dominates these data due to the quark charge 
factor, the positive Sivers amplitudes for u + (and h + , 
which is dominated by the pion yield), suggest a large 
and negative Sivers function for up quarks. The van- 
ishing amplitudes for tt~~ (h~) require cancellation ef- 
fects, e.g. from a d quark Sivers distribution opposite in 
sign to the u quark Sivers distribution. These cancella- 
tion effects between Sivers distributions for up and down 
quarks are supported by the vanishing Sivers amplitudes 
extracte d from deute r on da t a by the COMPASS C ollab- 



oration (|Ageev et all l2007t IXlekseev et all . l2009al) . An 



interesting facet of the HERMES data is the magnitude 
of the K + amplitude, which is near ly twice as large as 
that of 7r + ( Airapetian et al. . l2009bh . Again, on the ba- 
sis of u quark dominance, one might naively expect that 
the 7r + and K + amplitudes should be similar. Their dif- 
ference in size may thus point to a significant role of 
other quark flavors, e.g. sea quarks. A sizable Sivers am- 
pli tude for 7r + was also recently reported by JLab Hall 
A (jQian et all 120111) for measurements with a 3 He (neu- 
tron) target. In that data a negative Sivers amplitude for 
7r + was found which independently supports a d-quark 
Sivers distribution opposite in sign to the it-quark one. 

The Boer-Mulders distribution (Boer and Mulders . 



19981 ) describes the correlation between transversely po- 
larized quarks in an unpolarized nucleon and the quarks' 
transverse momentum, s q ■ (kt x p), where s q denotes 
the spin of the quark and might hence yield unexpected 
spin effects even in an unpolarized nucleon. It is similar 
to the Sivers distribution in that it is T-odd. However, 
it is also chiral-odd and hence must be probed in con- 
junction with a second chiral-odd function. For Drell- 
Yan production, the second function is a Boer-Mulders 
distribution in the second incident hadron. For SIDIS, 



the Collins fragmentation function described below is in- 
volved. Like the T-odd Sivers distribution, the Boer- 
Mulders distribution is also expected to change sign be- 
tween Drell-Yan production and SIDIS. Future experi- 
mental effort is planned to test this QCD prediction. 

Azimuthal distributions sensitive to the Boer-Mulders 
distribution were original l y mea s ured in Drell-Yan ex- 



periments (IConwav et all. 1989; 



Guanziroli et al. 



19881: 



Zhu et al. 



2007 



Falciano et al.. . 
2009). 



19861: 



The 



SeaQuest fixed-target Drell-Yan exp eriment currently 

7j) expects to be sen- 



2007 



underway at Fermilab ([Reimer 
sitive to the Boer-Mulders distribution at high x. In 
SIDIS the distinctive pattern of Boer-Mulders modula- 
tions for oppositely charged pions and f or pions and kaons 



was re cently reported by HERMES (|Airapetian et al. 
2012a ). The amplitudes for kaons are larger in mag- 
nitude than the amplitudes for pions. The amplitudes 
for the negative pions have the opposite sign to the 
amplitudes for negative kaons. This hints at a signif- 
icant contribution from sea quarks, in particular from 
strange quarks. Measurements of the Boer-Mulders am- 
plitudes were a lso repor t ed by COMPASS for uniden- 



tified hadrons (Sbrizzai|. | 201 lh and by CLAS for pi- 
ons ( Osipenko et all . 2009^ The interpretation of the 
SIDIS amplitudes for the Boer-Mulders distribution, is, 
however, complicated b y con t ributi ons from the twist- 
19781 Il989h 



4 Cahn effect (|Cahnl , 
estimated to be sizable 



which have been 
COMPASS kinemat- 



ev en at 

ics ( Anselmino et al. . 2007rJ ). The Cahn effect accounts 



for the parton intrinsic transverse momenta in the target 
nucleon and the fact that produced hadrons might ac- 
quire transverse momenta during the fragmentation pro- 
cess. Theoretical estimates of the Boer-Mulders effect 
are still plagued by large uncertainties, mainly related to 
the insufficient knowledge of the transverse-momentum 
dependence of the unpolarized distribution fi(x,kt) and 
fragmentation function D(z,pt). 



2. The Collins TMD fragmentation function 



COMPASS 2010 proton data 
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FIG. 17 Sivers amplitudes for unidentified ch arged hadrons 
measu red by COMPASS with a proton target ( Ador ph et all 
l2012bh . The hadron yield is dominated by pions. 



The Collins TMD fragmentation function describes a 
spin-momentum correlation in the hadronization process, 
s q • (k q x p t ), with a hadron produced in fragmentation 
having some transverse momentum p t with respect to the 
momentum direction k of a transversely polarized frag- 
ment ing quark with spin s q (ICollind . H99i ICollins et ~al\ . 
19941 ). The Collins fragmentation function has been 



investigated in semi-inclusive lepton-nucleon scattering 
and e + e~ annihilation. The magnitude of the effect is ap- 
proximately 5-10%, like that found for the Sivers asym- 
metries. 

For e + e~ annihilation the chiral-odd Collins fragmen- 
tation function enters with a second Collins function in 
the opposing jet. The Collins function has been measured 
to be non-zero for the production of charged pions in 
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FIG. 18 Collins asymmetry for the double ratios of like-sign 
(L), u nlike-sign (U) an d any charged (C) pion pairs from 
Belle (|Seidl et all l200St ). A UL and A uc are sensitive to dif- 
ferent combinations of the favored and unfavored Collins frag- 
mentation functions. 



e + e annihilation at Belle ( Abe et all , l2005HSeidl et al . 



2008), as show n in Fig.[T8l an d in recent preliminary data 



from BABAR (|Garzial 12012ft . 

In SIDIS the second chiral-odd function is the transver 
sity distribution introduced in Section II and dis- 
cussed further below (or the Boer-Mulders distribu - 
tion). The HERMES (fAirapetian et al.1. l2005bl |2010al) 



COMPASS lA dolph e t all 1201 2al lAgeev et 



lAlekseev et all 2009aL l2010bT iPesarol 1 2 1 lh and JLab 



2007 



Hall A (|Qian et all \20lW experiments have performed 
SIDIS measurements of the Collins effect. The measure- 
ments for a proton target are shown in Figs. [19] and [20] 
for HERMES and COMPASS respectively. (Note that 
COMPASS uses a definition of the Collins angle which 
results in Collins ampli tudes with oppo s ite sig n to the 
"Trento convention" of iBacchetta et all ([2004T ) used by 
HERMES, JLab and commonly in theoretical papers). 
There is excellent agreement between the measurements 
in similar kinematics. One finds the striking observation 
that the Collins amplitude for 7r~ is of similar size to tt + 
production but comes with opposite sign. This hints at 
an unfavored Collins function of similar size and opposite 
sign than the favored one, a situation very different from 
that observed with unpolarized fragmentation functions. 



3. Probing transversity 

The transversity distribution introduced in Section II 
describes the transverse polarization of quarks within a 
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FIG. 19 Collins amplitudes for charged pions measured by 
HERMES with a proton target (|Airapetian et all l2010al ) . 



COMPASS 2010 proton data 
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FIG. 20 Collins amplitudes for unidentified ch arged hadrons 
measu red by COMPASS with a proton target l|Adolph et all 
l2012bh . The hadron yield is dominated by pions. Note that a 
different definition of the Collins angle results in amplitudes 
with the opposite sign compared to other measurements. 



transversely polarized nucleon. Along with the unpolar- 
ized and helicity distributions, it survives integration over 
partonic transverse momentum and is thus a collinear 
distribution. 

The first moment of the transversity distribution is 
proportional to the nucleon's C-odd tensor charge, viz. 
Sq = J Q dxh\ (x) with 

(p, s\ qia^j 5 q \p, s) = (l/M)(s fJ ,p u - s„p^)Sq. (37) 

The difference between the transversity and helicity spin 
distributions reflects the relativistic character of quark 
motion in the nucleon. In Bag models this effect is mani- 
fest as follows. The lower component of the Dirac spinor 
enters the relativistic spin depolarization factor with the 
opposite sign to Aq because of the ex tra factor of 7 M 
in the tensor charge (|Jaffe and Jil . 1199211 . The relativis- 
tic Bag depolarization factor mentioned in Section VI 
becomes 0.83 for transversity in contrast to 0.65 for he- 
licity and the nucleon's axial- charges. In leading order 
QCD the transversity distributions are bound by Soffcr's 



inequa l ity \h q Ax,Q 2 )\ < \[{q + q}(x,Q 2 ) + Aq(x,Q 2 )], 



Sofferl ( 19951) . QCD m otivated fits to t r ansver sity ob- 
servables are reported in lAnselmino et all (|2009bh . which 
also reviews the comparison to model predictions. 

Transversity is measured through the Collins effect and 
also in in dihadron production, where the chiral-odd part- 
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ner of /if i s given by the d i hadron fragmentation func- 
tion H? q (Bacchetta et al. . 2011 ; Bianconi et al. . 2000t 



Collins et all 1994 ). This describes how the transverse 



spin of the fragmenting quark is transferred to the rel- 
ative orbital angular momentum of the hadron pair. 
Consequently, this mechanism does not require trans- 
verse momentum of the produced hadron pair. Standard 
collinear factorization applies allowing one to study the 
transversity distribution without having to worry about 
solving convolution integrals over transverse momentum 
or issues of TMD factorization and evolution. 

Pioneering measurements of two-pion produc- 
tion in polarized semi-inclusive DIS by HER- 



MES (lAirapetian et all. l2008al) and COM- 



PASS (jAdolph et all l2012dh reveal a sizable effect 



and have al ready been employed for an extraction of 



transversity (|Courtov et al 



2012a). First measurements 



of azimuthal correlations of two pion pairs in back-to- 
back jets in e + e~ annihilation related to the dihadron 
fragm entation funct i on ha ve just become available from 



Belle ( Vos sen et all 12011 ) and a first extraction of the 



dihadron frag mentation function f rom these data was 
performed in lCourtov et al. (|2012bl) . 



4. Current status and recent progress with TMD distributions 

There has been considerable progress in the under- 
standing of intrinsic transverse momentum and spin- 
momentum correlations in QCD over the past decade, 
motivated by the theore tical breakthroughs reg a rding T- 
odd T MD distributions (|Brodskv et aZ.I . l2002aiibl : ICollins . 



2002]) and by a vast program of theoretical and experi- 
mental activity. 



A recent monograph, ICollinsI (|2011l ) gives definitions 
of TMD distributions which allow QCD evolution to be 
applied rigorously for the first time with separately iden- 
tifiable TMD distributions and fragmentation functions. 
Building upon this progress, the evolution of previously 
unevolved models and fits has now been published for 
unpol arized TMD distributions and fragmentation func- 
tions J Avbat and Roeersl 2011bl ) and the Sivers distribu- 
tion (jAvbat et all 120121 ). QCD evolution i s just starting 



to be applied to phenomenological studies (jAvbat et al. 



201l[ ). which will be a major step forward in interpreting 
and comparing results from different experiments. The 
new definitions of TMD distributions also recently made 
possible a determination of the hard parts for SIDIS and 
Drell-Y an at next-to-leading order ( Avbat and Roge"rsl 
2011al ). which should lead to improved phenomenology. 

Much effort has been dedicated to phenomenolog- 
ical extractions of TMDs and parameterizations of 
the Sivers distribution from SIDIS data, see e.g. 



Anselmino et all (|2009al ) and lAnselmino et al. 



(201111 . 

with QCD evolution now starting to be consid- 
ered (jAnselmino et all I2012T) . One parameterization of 



the Sivers function includes both semi-inclusive deep 



inela stic and proton-proton data (jKang and Prokudin . 
2012). modulo issues re l ated to factorization break- 
ing (IRoeers and Muldersl . l2010l) discussed below. Fits 
to the Collins TMD fragmentation function have been 
per formed using both e + e~ and SIDIS data as in- 



put (jAnselmino et all l2007al lEfremov et all 120061) . The 



Boer-Mulders distribution has been extracted based on 



Drell-Yan (Lu and Schmidt. 201 



Zhang et all 120081 ) as 



2010 



well as SIDIS data (jBarone et all 

These first phenomenological fits to TMD observables 
have been performed using a simple Gaussian ansatz 
for the transverse momentum dependence of quarks in 
the nucleon and fragmentation functions. For exam- 
ple, the Sivers function in Eq.(36) was parametrized 
in the fits by the product of the unpolarized distribu- 
tion f q /pt(x,kt) with an x dependent factor and an x- 

independent Gaussian ~ jf^ e ~ kt ^ Ml containing all the 
kt dependence. While the Gaussian ansatz is unstable 
with respect to QCD evolution with increasing Q 2 , the 
method docs provide a reasonable fit to present data with 
values (fc t 2 ) = 0.25GeV 2 and (p 2 ) = 0.20GeV 2 taken 
from fits to the Calm effect in unpolarized scattering 
(jAnselmino et all |2005). A longer term goal for TMD 



experiments is to observe deviation from Gaussian behav- 
ior for transverse momentum dependence. With exten- 
sive unpolarized Drell-Yan and weak boson production 
data available over scales from ~ 4 GeV 2 to M§, new fits 
of unpolarized TMD distributions are quite promising as 
a means to test the Q 2 evolution of TMD distributions as 
well as to learn more about the shape of the distributions 
in kt. 

Lattice calculations of the Sivers and Boer-Mulders 



distributions have b e en performed ( Gockeler et all 12007 



Hagler et all 120091 : iMusch et all bfllj). There have 



also been efforts in recent years to i mplement TMDs 



in Monte Carlo eve nt generators (jBianconil . 12011 



Hautmann et all I2012T ). Models can provide helpful in- 



sight into TMD distributions, and a wealth of different 
model calcul ations have be e n exp lo red and pu b lished . 



We refer to lAvakian et all (l2009h. iBacchettal lEtt 



jitzerl 



Lorce and Pasauinil (|2011f ) , and iPasauini and Schweitzer! 
(|201lj ) for recent discussion of models related to TMD 
distributions, including attempts to address the relation- 
ship between TMD distributions and orbital angular mo- 
mentum in the nucleon 



5. Proton-proton asymmetries and TMD-factorization breaking 

Despite the fact that the large transverse single spin 
asymmetries observed in hadronic scattering originally 
inspired the development of TMDs, inclusive hadron pro- 
duction in p + p scattering cannot be cleanly separated 
into Sivers, Collins, or othe r contributions as is possi- 
ble in SIDIS. In recent work Rog ers and Mulderd (j2010l) 
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ANL 

^s=4.9 GeV 



BNL 

^s=6.6 GeV 



Fermilab 
ts=19.4 GeV 



RHIC u 
1^=62.4 GeV 



FIG. 21 The transverse single spin asymmetry in forward tt^ production as me asured in po l arized proton-proton c o llision s 
across a range of c enter- of-mass energ i es. F rom left to right, the data are from iKlem et al\ (1 19761 ) . lAllgower et al\ (|2002), 
lAdams et oZT (|l99lT) . and lArsene et al\ i|200Sl ). 



argue that the TMD framework is not valid in the case 
of hadroproduction of hadrons, as factorization does not 
hold. While the short-distance (perturbative) compo- 
nents are still believed to factorize from the long-distance 
(non-perturbative) ones, the long-distance components 
become entangled and no longer factorize from one an- 
other into independent TMD distributions and/or frag- 
mentation functions. What is particularly interesting is 
that the factorization breaking effects are relevant in pre- 
cisely the kinematic regime where a parton description is 
generally expected to apply. It will be exciting to see this 
tested experimentally in the upcoming years, exploring 
long-distance quantum entanglement effects in QCD. In 
the longer-term future, it may be possible to develop well- 
defined functions within the framework of pQCD which 
describe the correlations between the partons in the in- 
coming and/or outgoing hadrons. 

In the meantime, single spin asymmetries for forward 
meson production in p + p collisions have been shown 
to remain larg e across a very wide range of center-of- 



mass energies (lAbeley et al. . I2008ai l Adams e talx . 11991 



Il996t lAllgower et all l2002UArsene et aZj,l2008l) an d up 
to the highest measured pr of ~5 GeV ( Kosterl . 2012 ) . As 



shown in Fig. 1211 the transverse single spin asymmetries 
in charged pion production as a function of Feynman-ir 
are remarkably similar from y/s =4.9 GeV all the way up 
to 62.4 GeV measured by the BRAHMS experiment at 
RHIC. 

At higher energies and in particular at px values large 
enough to serve as a hard scale, one can try to inter- 
pret these phenomena utilizing the tools of pQCD. With 
no explicitly measured scale sensitive to the partonic 
transverse momentum in inclusive single spin asymme- 
tries, a more appropriate framework than TMD distri- 
butions in which to interp ret the asymmetries may be 
a col l inear, twist-3 picture dEfremov and Tervaev . 1982 , 



19851 : iQiu and Stermanl . I1999I ). A relationship between 



the TMD and the collinear, twist-3 frameworks was laid 



out in lJi et all (|2006j ). 

Surprises continue to emerge from these kinds of mea- 
surements, with large asymmetries f or negative kaons as 



well as antiprotons from BRAHMS (jArsene et all 120081 ) 



suggesting that the pion asymmetries are not a valence 
quark e ffect as previously believe d, and a recent hint from 
STAR (jAdamczvk et all l2012q) that the asymmetry for 
eta mesons may be larger than that of neutral pions. 



VIII. FUTURE PROJECTS 

A new program of dedicated experiments is planned to 
investigate key open questions in QCD spin physics. We 
briefly outline these experiments and their prime physics 
objectives. 

Since May 2012 CEBAF is undergoing a major up- 
grade that will bring the maximum available energy of 
the electron beam to 12 GeV. The experimental equip- 
ment in all three halls will be upgraded (Hall A and C) 
or completely renewed (Hall B), in order to better match 
the increased energy and luminosity. A new experimental 



Hall D is being built. Commissioning of the new acceler- 
ator and of the experimental halls is expected for 2014. 
The future physics program focuses on dedicated studies 
of large x phenomena, hard exclusive reactions and TMD 
effects i n kinematics where valence quarks dominate the 



physics (jDudek et al} , |2fll2). 



At CERN a proposal by the COMPASS Collaboration 
(|Gautheron et all l20irj ) to study TMDs and GPDs in 
the period 2014-2017 has been approved. The COM- 
PASS data will provide a link between the kinematic 
domains of HERA on one hand and of HERMES and 
JLAB on the other. The program will start with the first 
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ever polarized Drell-Yan experiment using a transversely 
polarized ammonia (proton) target and a negative pion 
beam. Due to the underlying annihilation of the anti up- 
quark from the pion and the target up quark, the process 
is dominated by the up quark distribution in the valence 
region. An important goal is to check the QCD prediction 
of a sign change in the naive T-odd TMDs with respect 
to the DIS case. A study of GPDs in DVCS and hard 
exclusive meson production with a polarized muon beam 
will follow in 2015 using a liquid hydrogen target, a dedi- 
cated target recoil detector and an additional large-angle 
electromagnetic calorimeter. An important measurement 
is the beam charge-and-spin asymmetry, which uses the 
property of the muon beam that polarization changes 
sign when going from positive to negative muons. A first 
result on the correlation of transverse size and longitudi- 
nal momentum fraction might already be expected from 
a 2012 pilot run. In parallel semi-inclusive DIS data will 
be taken on the pure hydrogen target. 

There are proposals to create a polarized fixed-target 
Drell-Yan program at Fermilab following the SeaQuest 
experiment, scheduled to complete data taking in 2014. 
R&D has begun for a suitable polarized target, and a for- 
mal proposal to polarize the 120 GeV proton beam in the 
Main Injector has been sub mitted to Fermilab manage- 



ment (jCourant et all 120111) . One of the primary physics 



motivations for such a program would be to explore in 
detail the QCD spin-momentum correlations described 
by TMD distributions such as the Sivers distribution, in 
particular the role of color flow in Drell-Yan versus semi- 
inclusive DIS interactions. 

A variety of possibilities for the medium-term future of 
RHIC is currently under discussio n. R&D is ongoing fo r 



a polarized 3 He source for RHIC (|Zelenski et aZ.I . [20081 . 



which would allow the neutron spin structure to be stud- 
ied in collider kinematics for the first time. There are 
also proposals to significa ntly extend the dete ctor capa- 
bilities at RHIC; see e.g. lAidala et al\ <|2012h . Of par- 
ticular interest to nucleon structure studies are potential 
upgraded forward spectrometers capable of reconstruct- 
ing jets, with hadronic particle identification and Drell- 
Yan measurement capabilities up to pseudorapidities of 
~ 4. The ability to perform full jet reconstruction in the 
forward region where large transverse single spin asym- 
metries are observed, and in addition to measure and 
identify hadrons within the jet, would allow separation of 
effects due to distribution versus fragmentation functions 
and shed great light on the origin of these significant spin- 
momentum correlations. An integrated design process for 
new detectors is underway such that they would be able 
to take full advantage of electron-proton and electron- 
ion collisions in the longer-term future should an electron 
beam be added to RHIC. 

Ideas for future polarization measurements are also dis- 
cussed and investigated in more detail at FAIR (Ger- 
many), J- PARC (Japan) and NIC A (Russia). 



A possible Electron-Ion Collider (EIC) is being dis- 
cussed in connection with the future of RHIC and JLab. 
The goal is to achieve highly polarized (greater than 70%) 
electron and light-nucleus beams with center-of-mass en- 
ergies ranging from about 20-150 GeV at maximum col- 
lision luminosities typically ~ 10 34 cm~ 2 s _1 . Signifi- 
cant R&D is ongoing to realize the technical challenges 
for reaching this luminosity frontier for colliders and for 
achieving and maintaining polarization of light nuclei (D 
and 3 He) in a storage ring. An EIC with the above per- 
formance would offer unique access to the small- a; re- 
gion where gluons dominate as well as to the intermedi- 
ate and high x regions at unprecedented high Q 2 . One 
could then st udy gluon polarization d own to x values of 



about 10 (jAschenauer et all 120121 ). The high lumi 



nosity would enable us to measure and map GPDs over 
a broad range of the kinematic variables and study the 
QCD evolution of the DVCS process plus TMD distribu- 
tions in kinematics where sea/glue effects are expected 
to be important. In addition to being the first ep col- 
lider exploring the structure of polarized protons, an EIC 
would also be the first electron-nucleus collider allowing 
precision studies of the gluon and sea quark structure 
of nuclei. Unpolarized ion beams from deuterium to the 
heaviest nuclei - uranium or lead - would also be ac- 
celerated. Knowledge about the spatial distribution of 
quarks and gluons in nuclei is needed for example in the 
interpretation of heavy-ion collision data and the search 
for quark-gluon plasma. A comp rehensive review o f EIC 
physics opportunities is given in lBoer et al. ( 2011 ). 



IX. CONCLUSIONS AND OUTLOOK 

The challenge to understand the internal spin structure 
of the nucleon has inspired a global program of enormous 
experimental and theoretical work in QCD during the 
last 25 years. 

For longitudinal spin structure, there is a good conver- 
gence of spin measurements from CERN, DESY, JLab, 
RHIC and SLAC taking into account the Q 2 dependence 
of the data and kinematics of the different experiments. 
There is also good convergence of theoretical understand- 
ing with the data, including QCD inspired models of 
the nucleon and lattice calculations with disconnected 
diagrams included. Semi-inclusive measurements in po- 
larized lepton-nucleon and proton-proton collisions have 
yielded much information about the size of the separate 
valence, sea and gluon spin contributions to the nucleon's 
spin. The small value of the nucleon's flavor-singlet axial- 
charge, about 0.35, extracted from polarized deep inelas- 
tic scattering seems to be a valence quark effect. No 
significant sea-quark polarization is observed in semi- 
inclusive deep inelastic scattering experiments; the sum 
of valence spin contributions is in close agreement with 
the measured total spin contribution <7^| p dis- While 
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gluon polarization Ag at the scale of the experiments may 
be as much as 50% of the nucleon's spin at the scale of 
the experiments, the QCD anomaly correction — 3§^A<7 

is too small to resolve the difference between g^| P Dis 
and the early quark models predictions, about 0.6. Prime 
theory candidates to explain the small "quark spin con- 
tent" include transfer of valence quark spin to quark or- 
bital angular momentum through the pion cloud and a 
possible topological effect whereby some fraction of the 
valence quarks' "spin" resides at Bjorken x = 0, where 
it is missed by polarized deep inelastic scattering exper- 
iments. The proton spin puzzle seems to be telling us 
about the interplay of valence quarks with chiral dynam- 
ics and the complex vacuum structure of QCD. Ongoing 
and planned experimental activity will improve the pre- 
cision on the size of gluon and strangeness polarization 
in the nucleon. 

Finite orbital angular momentum of the valence quarks 
is expected, induced also by confinement which intro- 
duces a transverse scale in the physics. Quark or- 
bital angular momentum through spin-orbit coupling is 
a prime candidate to explain the large transverse sin- 
gle spin asymmetries observed in proton-proton collisions 
and lepton-nucleon scattering. The desire to understand 
and measure QCD orbital angular momentum effects in 
the nucleon has spawned a new program to explore and 
map the three-dimensional structure of the nucleon - 
both in spatial co-ordinates (generalized parton distri- 
butions) and transverse momentum dependence. 

Studies of transverse nucleon structure will drive the 
experimental program in the near future, with dedicated 
running or approved programs at COMPASS, the 12 GcV 
upgrade of JLab, FNAL and RHIC. These experiments 
will test our understanding of initial and final state inter- 
actions in QCD (through comparison of Sivers and Boer- 
Mulders observables in Drell-Yan and semi-inclusive deep 
inelastic scattering). Precise measurements of GPDs and 
TMDs will test QCD evolution in a regime where trans- 
verse structure becomes important. The aim for precise 
information about quark (and gluon) total and orbital 
angular momentum in the nucleon is also a driving force 
for much theoretical work. Highlights include models of 
transverse spin phenomena, lattice calculations and de- 
velopment of QCD fitting technology to extract GPDs 
and TMDs from the newly measurable observables. 
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